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Abstract 

In this paper, we explain how perturbative quantum field theory can be formu- 
lated in terms of (a version of) vertex algebras. Our starting point is the Wilson- 
Zimmermann operator product expansion (OPE). Following ideas of a previous pa- 
per [arXiv: 0802 .21981 . we consider a consistency (essentially associativity) condition 
satisfied by the coefficients in this expansion. We observe that the information in the 
OPE coefficients can be repackaged straightforwardly into "vertex operators" and that 
the consistency condition then has essentially the same form as the key condition in 
the theory of vertex algebras. We develop a general theory of perturbations of the alge- 
bras that we encounter, similar in nature to the Hochschild cohomology describing the 
deformation theory of ordinary algebras. The main part of the paper is devoted to the 
question how one can calculate the perturbations corresponding to a given interaction 
Lagrangian (such as Xip^) in practice, using the consistency condition and the corre- 
sponding non-linear field equation. We derive graphical rules, which display the vertex 
operators (i.e., OPE coefficients) in terms of certain multiple series of hypergeometric 
type. 

1 Introduction 

It is becoming increasingly clear that Quantum Field Theory (QFT) is not only a very suc- 
cessful theoretical formalism to describe a variety of physical situations such as elementary 
particle collisions or critical phenomena, but also a highly interesting and complex math- 
ematical structure. This manifests itself not least in the range of different mathematical 
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disciplines that come to play an important role in this theory, such as functional analysis, 
combinatorics, complex analysis, algebra, geometry, measure theory etc. 

Not surprisingly, there exist correspondingly many approaches to QFT, such as the path- 
integral approach, the approach via diagrammatic/perturbative expansions, algebraic ap- 
proaches, stochastic quantization, axiomatic approaches, and more. These approaches are, 
at some level, known to be equivalent, but each one of them seems to have different weak- 
nesses and strengths. In a recent paper [H] (see also [15] for a related proposal), one of us 
proposed a new approach of an algebraic nature that is based on elevating the Kadanoff- 
Wilson-Zimmermann operator product expansion (OPE) [291 SO] to the status of a funda- 
mental structure. The operator product expansion is traditionally viewed as a property of 
the Schwinger functions of the composite fields of the theory. It can be stated as saying that 

i c i 

Here, x e is a point in Euclidean space not equal to 0, the yi are distinct points with 
\yi\ > \x\, and a,b,c etc. are indices labeling the composite fields of the theory, which 
in turn are monomials in the basic field and its derivatives. The OPE-coefficients C^f, are 
distributions that are independent of the "spectator fields" OdiiVi), whereas the brackets 
indicate the Schwinger functions of the (Euclidean) QFT under consideration. They may be 
constructed e.g. by giving precise mathematical sense to a functional integral of the type 

(YlOaAx^d) = I WOaM) exp[-5A(vP)] d/iA(^). (1.2) 

either in perturbation theory [21], [22], or non-perturbatively, see e.g. [3 [TBI [25], [23]. In the 
above formula, is a Gaussian measure associated with the free part of a classical 

action, and is the full action including all counterterms, depending on the cutoff scale A. 
The OPE coefficients characterize the short-distance properties (|x| 0) of the Schwinger 
functions. The formula is normally understood as an asymptotic expansion in this limit, 
with coefficients becoming more and more smooth in their argument as the dimension 
of the operator indicated by Oc increases. 

The proposal of [lU [16] was to view the OPE coefficients not as a secondary structure, 
but as the defining structure of the theory — some kind of "structure functions" of a suitable 
algebraic structure governing the field theory. In order to display this algebraic structure, 
we looked at Schwinger functions with a triple insertion Oa{x)Oh{y)Oc{Gi), assuming that 
< |a; — |/| < |?/| < It is then plausible that one obtains equivalent results when 
performing successive OPE's in the following two alternative ways indicated by the different 
positioning of the parenthesis in (Ca(a;)0b(|/))0c(O) = Oa{x){Ob{y)Oc{Q)). For example, 
the left side means that we perform an OPE of the fields Ofy{y) with the field Oc{0), and a 
subsequent OPE of the result with Oa{x). The equality sign between both sides means that 
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the results should be the same when inserted in a Schwinger function as above. In terms of 
the OPE-coefficients this means that we expect the relation 

Y^Ciix - yWAy) = Y.Ct{y)CUx) . (1.3) 

d d 

We called this the "consistency" or "associativity" condition in P^j- We observe that this 
relation becomes much more transparent if we introduce the following notation: first, we 
introduce an abstract infinite dimensional complex vector space, V , whose basis elements 
a & V are in one-to-one correspondence with the composite fields Oa of the theory. Thus, 
by a slight abuse of notation, we use the labels a, 6, c both as labels of the fields, and 
simultaneously for the basis elements of V . We also define the corresponding (algebraic) dual 
vector space V* . We then define a "vertex operator" Y {a^x) : V ^ V a.s the endomorphism 
of V whose matrix elements are given by 

(c|F(a,x)|6):=C:,(x), (1.4) 

where we use the usual physicist's notation for basis vectors of V as "kets" \a) and the 
corresponding dual basis vectors of V* as "bras" With this notation, the consistency 
condition reads simply 

F(a, x)Y{h, y) = YiY{a, x - y)b, y) , (1.5) 

which holds again when |x| > \y\ > \x — y\ > 0. In this form, our condition is almost 
identical in appearance to (one form of) the key relation in the theory of vertex operator 
algebras [IHl El [9l |TT] , hence the name "vertex operator" for Y. We explain the relationship 
of our approach to vertex algebras in somewhat more detail in sec. [6] of this paper, but 
we emphasize from the outset that our setup is intended to be much more general than 
that usually encountered in this theory; for example, we do not assume that the theory is 
conformally invariant, and we admit general dimensions D, even though we expect non- 
trivial examples only in D < 4. As a consequence, the contents of our consistency condition 
and that in vertex algebras are actually rather different despite the similar appearance. An 
alternative treatment of the free field vertex algebra in D dimensions can already be found 
in ^ and [1]. 

The main purpose of the paper is to outline how the standard constructions in QFT, 
such as perturbative expansion, BRST cohomology etc. can be formulated in terms of the 
vertex operators Y, and what new viewpoint one thereby gets for these constructions. As 
we will see, if we combine the consistency condition with a perturbative expansion, and with 
the non-linear field equations corresponding to the classical action in the path integral 
above, then we obtain in a new scheme for doing perturbative calculations. In more detail, 
the contents of the paper are as follows: 

• In sec. |2]we begin by explaining the general setup. In particular, we explain in abstract 
terms how ordinary perturbation theory fits into the framework of vertex algebras. 
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As we show translating constructions of [E] into the framework of vertex algebras, 
one can define a cohomology of "Hochschild-type" associated with a D-dimensional 
vertex algebra (in our sense), and perturbations may be characterized in terms of this 
cohomology. We also explain how the framework may be generalized to also include 
the BRST-construction for gauge theories. 

• In sec. El we explain how the vertex algebra of the D-dimensional free massless field 
can be obtained from the Schwinger functions, and we thereby obtain a concrete (albeit 
somewhat trivial) example for our setting. This section also serves as the starting point 
to the later sections. 

• In sec. m we then construct the deformations ( "perturbations" ) of the free field vertex 
algebra corresponding to an interaction term such as Xip^ in the classical Lagrangian. 
We show how one can obtain a closed form expression for the i-th order deformation 
of the vertex operators based on the repeated use of the consistency condition and the 
(non-linear) field equation [see eqs. (15.1141) and (IS.llip for the final form of these ex- 
pressions]. The case D = 2 is in many ways simpler than the case D > 2, so we discuss 
it separately. However, in both cases, the final result is expressed in terms of sums 
that are associated with certain tree graphs and associated loop graphs. The resulting 
sums and integrals involve Gegenbauer/Legendre functions and are reminiscent of the 
expressions one obtains using the "Gegenbauer polynomial x-space technique" of [6j 
for multi-loop Feynman diagrams. 

• In sec. O we discuss an alternative representation of the deformed vertex operators 
in terms of certain multiple infinite sums. These sums have a similarity with the 
hypergeometric series, and we believe that this opens up an interesting connection 
between our vertex algebras and a class of special functions. 

• In sec. [6] we discuss in some detail the similarities and differences between our notion of 
vertex algebra and notions that were previously given in the context of 2-dimensional 
CFT's. 

• Various results and definitions related to Legendre functions in D-dimensions needed 
in the main text are given in app. \M and [Bi 

This paper is about the mathematical structure of QFT, but it is by itself not com- 
pletely mathematically rigorous. We analyze the structure and exploit the consequences 
of the consistency condition in a mathematically rigorous way, but we do not justify 
the consistency condition itself (apart from free field theory) and rather treat it as a 
hypothesis. In this sense, our approach can be thought of as some sort of "bootstrap". 
In order to justify the consistency condition, one can e.g. start from the OPE of the 
Schwinger functions, and show that performing the expansion of a triple product of fields 
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lead to the same result, as we described. We will come back to this issue in another paper [T^ . 

Notations and conventions: x,y etc. denote points in M.^, with scalar product x ■ y = 
^^^nVtii and norm = |xp = x ■ x. We also use the notation x := x/r G S^~^ for the 
angular part if x 7^ 0. P(2, z/, D) denote the Legendre functions, see appendix El and ip = 
r'/r is the Psi-function. The surface area of the {D — l)-dimensional sphere is abbreviated 
an = yId/2) = V-D — 2. The natural numbers N include 0, N = {0, 1,2,...}. 

2 General setup 

In the previous section, we explained the basic idea that connects QFT as formulated in terms 
of Schwinger functions with vertex operators Y{a,x) satisfying a consistency condition. Let 
us isolate the properties that we expect these vertex operators to have, reflecting the general 
known properties of the Schwinger functions: 

1. First, we state formally that there is an identity operator in field theory. This is a 
distinguished element 1 & V . The OPE of 1 with any other field is trivial and the 
vertex operator associated with 1 is hence given by Y{l,x) = id, where id is the 
identity on V. 

2. We have already said that the OPE coefficients are in general expected to be distribu- 
tions on M^. In Euclidean space, the situation is actually better and the coefficients 
are analytic except for the origin. Thus, if we denote by ^(M^ \ {0}) the set of such 
analytic functions then we expect that Y{a, . ) G ^(M^ \ {0}) (g) End(l^). 

3. The Schwinger functions have an obvious invariance property under the Euclidean 
group SO{D) (or its covering group, if there are spinor fields in the theory). 
These invariance properties are inherited by the OPE coefficients, and we hence 
have corresponding invariance properties of the vertex operators: Y{a,g ■ x) = 
R{g)Y{R{g)~^a, x)R{g)~^ , g G SO{D), for some representation R of SO{D) on V. 
Note, however, that there is no similar covariance under translations, as the vertex 
operator in effect depend on the choice of the origin G as the reference point in 
the OPE. 

4. The key condition is the consistency condition, which we have already motivated and 
which we repeat: 

y (a, x)Yib, y) = F(r (a, x - y)b, y) , (2.6) 

for all a, 6 G y and all x,y E subject to < \x — y\ < \y\ < \x\. Note that there is 
an implicit statement about the convergence of an infinite sum made here: If we apply 
e.g. the right side of the equation to a vector c & V, then we effectively state that 
Y{b,y)c is in the domain of Y{a,x). When written in a standard basis of composite 
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fields, the vector d = Y{b,y)c typically is an infinite linear combination of such fields 
even if b, c are from a basis of composite fields. Thus, Y{a, x)d will also be an infinite 
sum when written in the basis of composite fields, and we require, in effect, that this 
sum converges. As we discuss in more detail in sec. [6], there is no such issue in the 
usual formulation of vertex algebras, where all sums are either formal, or finite. 

The Schwinger functions have other properties that one can readily translate into prop- 
erties of the vertex operators, such as (anti-) symmetry properties under the exchange of 
fields, scaling properties, hermiticity properties, etc. To keep the discussion transparent at 
this point, we will not go into this here, but state some of these conditions as we go along. 

In renormalized perturbation theory, one is naturally led to consider field-redefinitions, 
for example when considering changes in the renormalization scheme. A field redefinition 
is simply a linear transformation which maps a quantum field to a linear combination of 
quantum fields, and thereby gives a corresponding transformation of the Schwinger functions. 
Such field redefinitions change in an evident way also the OPE coefficients of the theory. In 
terms of the vertex operators, a field redefinition is simply an invertible complex linear map 
Z : V ^ V which commutes with R{g) for all g G SO{D), and which satisfies Zl = 1. 
The OPE coefficients associated with the transformed fields give rise to a correspondingly 
transformed vertex operator, given by 

Y'{a,x) = ZY{x,Za) Z-^ (2.7) 

This vertex operator satisfies the same properties as above, and will be considered "equiva- 
lent". 

An important construction in QFT is to consider perturbations of a given theory. The 
idea is usually to start from a free field theory with explicitly known Schwinger func- 
tions/OPE and to construct perturbations of these quantities, order by order in a "small" 
parameter A characterizing the size of the perturbation. In the path integral, one expands 
exp(— S'a) in terms of this parameter, and then performs the Gaussian integral for each term 
in this expansion. Such a procedure will also give rise to a corresponding expansion of the 
OPE coefficients, see e.g. |20J, and thereby to an expansion of the vertex operators of the 
form 

oo 

y(a,x) = ^AVi(a,a;). (2.8) 

1=0 

The zeroth order term Yo{a, x) of this — in general only formal — expansion corresponds to the 
free field theory around which one expands. The higher terms represent the perturbations. 
Because the full Y{a,x) must satisfy the properties 1) — 4) stated above for all A (or at least 
in the sense of formal perturbation series), the expansion coefficients Yi{a,x) must satisfy 
analogous properties. The most prominent one is the consistency condition, which at i-th 
order reads 

i i 

x)Y,_,{b, y) = Y, yAYi-A<^^ ^ - y)^, y) ■ (2.9) 

j=0 j=0 
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It is satisfied at order i = 0, because one can prove it for a free field tlieory. As already 
pointed out in [Hj, there is an interesting cohomological interpretation of eq. 02.91) . and we 
now briefly explain how this comes about. Let n > and define Q'^~^^{V) to be a space of 
linear maps 

...,Xn) -.V (g)---(g)V ^ End{V) (2.10) 

that are defined for certain configurations (xi, . . . ,Xn) G M"'^ of mutually different points, 
and where there are n tensor copies of V. More precisely, /„'s should be real analytic 
functions of Xi, . . . , x„ G M"^ in the open domain 

Dn = {{xi,...,Xn) e M"^ : ri,,_i < r^^i^i < ri„2,i < ••• < n^i] , (2.11) 

where Vij = \xi — Xj\, taking values in End(V^'^", End(V^)). Next, we define from Yq a linear 
operator b : f2"(V^) — >• f2"+^(V^) by the formula 

{bf n){xi, . . . ) •— ^o('^l7 2;i)/„(x2, . . . , Xn+l] a2, . . . , fln+l) 

n 

+ ^(-i)7„( Xi, . . . , Xi, . . . , Xn+l] 0,1, ... , lo(C'j, Xi — Xj_|-i)aj_|-i, . . . Ctn+l) 
i=l 

+ (-l)"+Vn(Xl,... ). (2.12) 

This linear operator is presumably not defined on all such /„, but only on a certain domain 
of definition, but we ignore this issue here for simplicity. We formally calculate [Hj using 
the consistency condition for Yq that 

W„) = 0, (2.13) 

i.e., 6 is a coboundary operator. The domains Dn are needed in order to be able to apply 
the consistency condition. We denote by 

H''{V,Yo) := {kerb: 1]" ^ n"+^}/{mnb : Q""'^ Q"} (2.14) 

the n-th cohomology ring of the chain complex (©„fi"(y),6). 

The connection to the consistency condition eq. (12.91) at i-th order in perturbation theory 
now arises as follows: At i = 1, the equation can be restated as simply saying that 

bYi = 0. (2.15) 

If Yi arises from Yq by merely a field redefinition (12.71) with Z = XliLi -^j-^N then this means 

Yi = bzi. (2.16) 

Consequently, the non-trivial first order perturbations may be viewed as elements in 
H'^{V,Yo), and hence such perturbations only exist if this space is non-trivial. Continu- 
ing in this way, the condition (12.91) at i-th order can be stated as saying that 

bYi = Wi , (2.17) 
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where Wi G fl^iV) is defined by the terms in eq. fll.3p with j 7^ 0, j 7^ i. It can be checked 
inductively [H] that bwi = 0, so [wi] defines a class in H^{V,Yq). Evidently, the i-th order 
perturbation Yi exists, if, in fact, Wi G ran 6, i.e. if it defines the zero class in H^{V, Yq). The 
freedom of choosing different solutions to eq. (12.171) corresponds to elements in H'^{V, Yq). In 
this way, we get an interpretation of perturbation theory in terms of the cohomology rings 
H-{V,Yo). 

There is no difficulty to extend this framework to more complicated situations where 
one has additional symmetries that one would like to preserve in the deformation process. 
Suppose e.g. that the free theory is a free U{1)^ gauge theory with associated ghosts. Then 
in the free field theory we have an associated free BRST-operator sq with the property Sq = 
and SQogh = (gh — id) osq, where gh is the "ghost number", i.e. a linear operator on V with 
integer spectrum that provides a corresponding grading. This acts on the composite fields 
of the theory, and hence, as a linear transformation Sq : V ^ V satisfying 7S0 + So7 = 0, 
where 7 : ^ — > ^ is a Z2-grading of V that corresponds to the Bose/Fermi character of the 
fields. The concrete form of sq in terms of the gauge field and ghost /auxiliary fields is given 
in many textbooks, see e.g. [31j- To explain our general scheme it is only important to know 
that the BRST-invariance of the Schwinger functions of the free gauge theory gives rise to 
a corresponding invariance of the OPE, namely, the OPE of BRST- invariant fields again 
contains only such fields. In terms of the vertex operators, this is expressed by the relation 

SoYo{a, x) + lo(7a, x)sq + lol^oa, x) = . (2.18) 

We would now like to ask whether it is possible to find a deformation of Yq as above to a F , 
and a corresponding deformation of sq to an s 

00 

s = 5^s,AS (2.19) 

such that = 0, i.e. 

i 

= ^s,Si_,-, (2.20) 
and such that eq. (I2.18P continues to hold for Y, s, that is, at i-th order, 

i 

SjYi^j{a, x) + Yj{'ja, x)si-j + Yj{si-ja, x) = . (2.21) 

j=0 

The problem of finding simultaneously the Si and Yi subject to eqs. (I2.2ip . (I2.20p and (12. 9p 
is again of a cohomological nature. To set up the cohomology ring in question, we extend 
the action of Sq from V to Q"'(y) by setting 

{Bfn){xi, . . . ,Xn;ai, . . . ,an) := sofn{xi, . . . ,Xn;ai, . . . ,an) (2.22) 

n 
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Then one verifies [13] that = using Sq = and that bB + Bb = using eq. (12.181) . i.e. B 
is another coboundary operator which is compatible with b. In the language of differential 
complexes, (6, B) gives rise to a double complex {(Bn,g^"''^ {V) , b, B), whose cohomology rings 
are denoted H^'^{V, Yq, sq). Here n is as above, and g corresponds to the grading of V (and 
correspondingly the spaces fl"'{V) = (Bgfi"''^iV)) by the ghost number gh. We have the 
following commutative diagram: 



B 



H'^%V, Yo, So) ^ H''\V, Yo, So) H'^'{V, Y, 



B 



0) So) 



B 



B 



B 



H''\V, Yo, So) H''\V, Yo, so) — H''\V, Yo, so 



B 



B 



As is standard in this situation, we may form the total complex whose cohomology rings are 
given by 

R'^{V,Yo\so) = H^'^{V,Yo,So), (2.23) 

n+g=m 

and whose total coboundary operator is + i.e. {B + b)'^ = 0. The direct sum in eq. (12.231) 
consists of the cohomology groups lined up on the diagonal running from if™'°(V, Yo, So) to 
if°'™(l^, Yq, So) in the above diagram. 

If we form = (sj, 1^, 0, 0, . . . ), then the combined condition for the first order perturbation 
ai arising from the first order consistency condition and the first order BRST-invariance can 
be written as 

{b + B)ai = 0, (2.24) 

while the condition for the first order perturbations to be due to a first order field redefinition 
(i = {zi, 0, 0, ... ) is written as 

ai = (5 + b)Ci . (2.25) 

Thus, first order perturbations of the BRST operator and vertex operators are given by a 
class [ai] G H'^{V, Yq\so). Likewise, the conditions for the z-th order perturbations a, can be 
written as 

(B + b)ai = Pi , (2.26) 

where (3i is calculable from the aj, j < i — 1, and where one can compute {B + b)j3i = 0. 
Thus, the potential i-th order obstruction is the class [Pi] G H^iV, So\Yo). The details of this 
analysis are completely analogous to [H]. 
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3 The free field OPE vertex algebra 

In this section, we illustrate our abstract framework for the OPE in a simple example. Our 
example is the free quantum field theory obeying the linear field equation 

Aip = 0. (3.27) 

The space V of fields in this theory may be taken to be the unital, free, commutative 
ring generated by the identity (/? and its derivatives. In other words, the elements of V 
are in one-to-one correspondence with monomials in . . - d^^ip, and d^,^ = 1, ...,D are 
the derivations that act as if they were ordinary partial derivatives. To implement the field 
equation, we simply set to zero any expressions containing a factor the form d^^^^d^j^^ . . . d^^^ip, 
i.e., monomials that would vanish if Lp was an actual field satisfying the field equation. 
Because monomials containing a trace of d^^ . . -d^^ip are set to zero, V is spanned by all 
trace-free monomials. Thus, if we denote by curly brackets t^^^,,,^^,y the trace-free part of a 
symmetric tensor, then a basis of V is given by 1, together with the set of monomials of the 
formn'9{/.i---<9^,}<^. 

It is convenient for latter purposes to choose a particular basis. For this, we consider the 
space of harmonic polynomials in D real variables homogeneous of degree Z, i.e. the set of all 
polynomials h{x) in D variables (x G MP) with complex coefficients satisfying h{tx) = t^h{x), 
and Ah{x) = 0. Some relevant facts about such polynomials are collected in appendix|Xl We 
denote by hi^m, m = 1, . . . , N{D, I) a basis of degree / harmonic polynomials. (The number 
of linearly independent polynomials of this kind, N{D, I), can be found in appendix lAl) We 
normalize this basis so that 

/ hi^rn{x)hi>^rn'{x)dQ{x) = 6i^i:6m,m' (3.28) 

where dfl is the standard integration element on the sphere, and x = x/\x\. A basis of V is 
then given by 1, together with the elements 

a = ll^={c;'hUd)^y'- , (3.29) 

l,m V^l,^- 

where a = {ai^rn € N | Z > 0, m = 1, N{1, D)} is identified at the same time with a vector 
in V and a multi-index of non-negative integers, only finitely many of which are non-zero, 
and ci is given in appendix O The Schwinger functions of the model are well-known. For 
2s factors of the basic field they are given by {D > 2) 

2s s 

Y[<p{Xi)^ = Yl\Xa{2j-l) - Xa{2j)\^~^ (3.30) 

i=l o-eSym{2s) j=l 



^With this normalization, the harmonic polynomials restricted to ^ are the _D— 1-dimensional spherical 
harmonics. 
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for non-coincident Xj G M^. For 2s + 1 factors of the basic field the Schwinger function is zero. 
Partial derivatives can be taken on both sides of the equation to get the Schwinger functions 
of all (yj-linear fields. The Schwinger functions of composite operators can be obtained in 
the same manner, writing a composite field as a product of yj-linear fields with the same 
argument x G M.^ and omitting permutations a with a;o-{2j-i) = Xa{2j) for some j G {1, . . . , s} 
in eq. (13.301) . For example, if none of the composite field has derivatives, we have 

n 

nv^""w)= j2 n i^v-x^i'-"". (3.31) 

v=l graphs G e=(tiui)eG 

Here, the sum is over all graphs with coordination numbers {pv}, and the product is over all 
edges e = (vw) of the graph. For composite fields with derivatives, one has a similar formula. 
From the Schwinger functions, one gets the OPE coefficients according to eq. (II. ip . and from 
the OPE coefficients, one gets the Yo{a,x) according to eq. (11.41) . In D = 2 dimensions, the 
propagators are replaced by In \xy — x^l- We now present the results of these calculations, 
leaving the details to appendix O The subscript "0" reminds us that we are dealing with a 
free field in this section. 

In order to present explicitly the lo's, it is convenient to view \^ as a "Fock-space," with 
ai^rn (see eq. (13.291) ) interpreted as the "occupation number" of the "mode" labeled by l,m, 
and with 1 playing the role of " Fock- vacuum" denoted |0) (vanishing occupation number). 
On this Fock-space, one can then define creation and annihilation operators hi^^, '■ ^ ~^ 
V, see appendix O They satisfy the standard commutation relations 

[\m,K,rn'] = ki'^m,m' id, [b+„, b+^,] = [b^,™, b;,_„,] = (3.32) 

where id is the identity operator on V". In this language, the basis elements of V are written 

as 

« = nM^|0)- (3-33) 

l,m V"'-'"' 

We now give the formula for YQ{ip,x) corresponding to the basic field. For D > 2, this is 
given by 

oo N{l,D) 



Yo{v,x)=Knr-^^-'^/^J2 E 



1=0 :mJl V^iDJ) 

V'+(^-^)/^/.,™(x) + r-^-^^'^y^hU^)h,J , (3.34) 



where Kd = \JD — 2, and the "frequency" co'(/, U) is given by 2Z + D — 2, see appendix O 
For a general a G of the form eq. (I3.29p . the vertex operators yo(a, x) :V are 

Yo{a,x) = : n r m/2 {c.^Krrmy^i^.AT''^ ■■ ■ (3.35) 

l,m ' 
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Here, double dots : ■ ■ ■ : mean "normal ordering", i.e., all creation operators are to the right 
of all annihilation operators. 

It can be checked explicitly that the Iq satisfy the properties 1) — 4). The representation 
R : SO{D) —>■ End(\^) in 3) is characterized by 

N{D,l) 

R{g) K„^R{gr' = '^L^m'igHm' ' (3-36) 

m'=l 

where in turn is the representatioij^ of SO{D) on the space of degree / harmonic poly- 
nomials, i.e. J2m''^'mm'i9)^im'i^) ~ mid ' ^) ■ Order to demonstrate the consistency 
condition 4), one must make use of the identities related to the harmonic polynomials given 
in appendix El We are not going to go through this lengthy calculation here. 



4 Perturbative OPE vertex algebras 

In the previous section we have described the (well-known) OPE of the free, massless bosonic 
field in D dimensions. The new viewpoint was to understand this as a vertex algebra. Of 
course, free field theories are only of limited interest, and we therefore now look at per- 
turbations of the free theory. Normally, perturbations are characterized by an interaction 
Lagrangian, which, together with appropriate counterterms, is inserted into the path inte- 
gral (11.21) in order to obtain the perturbation series of the Schwinger functions. To make this 
work, one first considers a regulated path integral, with a regulated interaction, S\, and one 
then removes the regulator A. This procedure is explained in many textbooks, see e.g. [31], 
and it leads also to the definition of the OPE coefficients, C^^, see e.g. [20] for a derivation 
using the Wilson- Wegner-Polchinski RG-fiow equations. 

As we have explained, in this paper, we want to pursue an approach wherein the OPE 
is elevated to the status of a fundamental relation, and we should therefore also have a 
method to calculate the perturbations of the OPE coefficients — or equivalently the F's — 
directly, without recourse to the Schwinger functions. In principle, we have outlined how 
this works in sec. [21 But for this we would need to understand more explicitly the — rather 
abstractly defined — cohomology rings H'^iV, Yq) and H^iV, Yq). We have so far been unable 
to do so. However, one would think that perturbations should also be describable in terms of 
an interaction Lagrangian, in the same way as in the usual approach, because one is talking 
about the same theory after all. 

We now explain one way how to proceed. Let us assume that the "bare interaction" 
(in the usual QFT parlance) is given by a polynomial XP{^p) = A ^ ^V^^- ■'■^ order to get 
a well-defined perturbative definition of the Schwinger functions, we also assume that the 

^In D — 3, I corresponds to the angular momentum quantum number and is the corresponding 
irreducible representation. For D > 3, the representation X>' is reducible. 
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interaction is renormalizable, i.e. deg(P)-^2^ < D. It is a well-known fact in standard 
perturbation theory that the Schwinger functions of the theory may then be defined so that 

([Aifix) - XP'{ip{x)) PM n 0,Xy^)) = . (4.37) 

i 

Here, Oa, Od, are arbitrary composite fields, and the arguments satisfy \yi\ > \x\ > 0. I The 
expression in brackets [. . . ] is of course the non-linear field equation. Let us now apply the 
OPE to the expression [Aip{x) — XP'{ip{x))]Oa{0) in the above Schwinger function. Then 
we get a relation between the OPE coefficients involving Aip and those involving P'{^p)- In 
terms of the vertex operators this relation is 

AY{ip,x) = XY{P'{ip),x) , (4.38) 

where we are now viewing P'{^p) as an element in V, the abstract vector space whose elements 
are in one-to-one correspondence with the composite fields. 

Next, we expand the vertex operators in a (formal) perturbation series in A as in eq. (12. 8p . 
and this immediately leads to the relation 

AF,((^,x) =l^_i(P'(v9),x). (4.39) 

The evident strategy is now to try to design an iterative scheme from this equation, by 
calculating the order i vertex operator on the left in terms of the lower order i — 1 vertex 
operator on the right, starting inductively with i = 1, as all vertex operators of order i — l = 
have been given in the previous section. However, such a procedure runs into the immediate 
difficulty that the right side involves the vertex operator associated with the composite field 
P'{(p), whereas the left side only gives the vertex operator of the basic field (p. We must 
therefore introduce a second induction loop which constructs, iteratively, the vertex operators 
of an arbitrary a & V from those of if. It is the essential strength of our approach that this 
is possible, using the consistency condition in perturbative form, see eq. (12.91) . 

Thus, suppose that we are given Yi{ip,x), and all other vertex operators up to order 
1, . . . , i — 1. How to find Yi{(p^, x), Yi{ip^, x), etc.? Consider points such that < \x — y\ < 
\y\ < \x\ and the following special case of the i-th order consistency condition: 

i i 

J2yJ{v,^)y^-Av,y) = Y.^^~J{yl{v,x - y)v,y) ■ (4.40) 

j=Q j=0 

Under the hypothesis that this condition indeed holds for the — yet to be constructed — terms 
in this equation that are not known at this point, we can draw the following conclusion. Let 
us look at the term with j = on the right side of the equation. Using the known form of 
the Yo's (from the free theory), we have Yq{lp, x — y)(f = \x — y]^'-'^^^-' 1 + v^^, plus terms that 



Otherwise, there would also be "contact terms" in the above equation. 
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are smooth in x — y and vanish for x = y. Using ^^(l, y) = for i > 0, we hence arrive at 
the relation 

i i 

Yi{ip\ y) = J2 x)y^~J{v, y)-Yl ^ - y)^' !/) + •••• (4-41) 

i=o j=i 

The dots represent terms that vanish in the hmit as |x — ?/| 0. The singular terms with the 
— signs may be thought of as some sort of "counterterms" , which cancel off the divergence 
of the first term on the right side in the limit. The key thing to note is now that all terms on 
the right side are known, by induction. Thus, taking the limit, we obtain the desired vertex 
operator Yi{ip'^,x), and, by iterating this procedure in an obvious way, all other i-th order 
vertex operators Yi{ip^, x), Yi^ip"^, x), . . . . 

In summary, our iterative scheme is set up in such a way that, at order i, we have to 
perform one inversion of the Laplace operator in eq. (14.391) to get Yi{(p,x), and then we 
subsequently have to construct all other vertex operators at order i via the consistency 
condition. Unfortunately, it is not evident from what we have said that the vertex operators 
constructed in this way will satisfy the consistency condition to all orders. Furthermore, 
there is certainly the freedom to add to the i-th order solution to the inhomogeneous Laplace 
equation Yi{ip, x) a solution to the homogeneous equation. It would seem that both issues are 
connected, and that one must impose the validity of the i-th order consistency condition in 
order to (partly) fix this ambiguity. One would then, furthermore, expect that the ambiguity 
is equivalent to the usual sort of renormalization ambiguity, which in our framework is given 
by eq. (12. 7p (with Z = ^ A*Zi) and a change in the coupling constant A by a (formal) 
diffeomorphism. 

We will not address these issues in the present paper but rather focus on the kind of 
mathematical expressions that one obtains following the above iterative scheme. We divide 
our discussion into several parts. First, in the next subsection we discuss how to take the 
inverse of the Laplace operator in a way that is suitable in our setting. Then in subsec. 14.21 
we will discuss how to organize the terms that appear in the iteration process by a graphical 
notation involving trees. In subsec. 14. 3[ we outline how the "counterterms" may be incor- 
porated into the graphical notation. In sec. [H we find another representation of the vertex 
operators in terms of infinite sums of hypergeometric type. 

4.1 Inverting the Laplace operator 

In our inductive scheme, we have to invert the Laplace operator at each step of the induction. 
According to the general setup, the vertex operators are functions Y{a, . ) G A^K^ \ {0}) (g> 
End(V^), so we should invert on this function space. However, in perturbation theory, the 
space of functions is actually much more restricted. At zeroth order, the vertex operators are 
in fact infinite sums of products of creation/anhihilation operators, harmonic polynomials 
hi^mix) and powers r^, where / gN, /cgZ, 0<m< N{D,l), see eqs. (I3.34p and (I3.35p . 
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Unfortunately, such functions are not stable under the inverse of the Laplace operator — we 
also get factors of Inr. Inverting again, we get factors of In^r, and so on. To incorporate 
the logarithms, we are thus led to introduce the following spaces of functions for i G N: 

£i = span {r'' W r hi^m{x) : A; G Z, j, / G N, j < i, 1 < m < A^(/, D)] . (4.42) 

Evidently, the union Uj£j is a filtered ring {£i£j C £i+j)- Any right inverse of the Laplace 
operator, G, takes us between £i £i+i for all i. Our inductive scheme will hence give us 
2-th order vertex operator in the space 

Yi{a, .)e£i® EndiV) . (4.43) 

In order to give an explicit formula for G we introduce a residue integral trick for computing 
right inverses of the Laplacian which we are going to use extensively. We first define 

27rifc 5 {k + 2 + S){k + D + S)-l{l + D- 2) ''"^"'^ ' ^ ^ ^ 

This defines G as an operator on £o —>■ £i. To extend the trick fl4.44p to all of £^ = ^j£j, 
assume that we are given f{x) G £j as a j-fold residue integral of the form 



,5,)r''hi,^{x) , (4.45) 



where F{5i, . . . ,5j) does not depend on x. It is always possible to represent / G £j as a 
linear combination of expressions of the form (14.451) . and in fact, this is the form we will 
encounter. We define the Laplace inverse of f{x) by 



X 



j+1 i+1 
{k + 2 + Y,5i){k + D + Y^5i)-l{l + D-2) 

i=l i=l 



-1 



r'^'K^{x).{A.m) 



We remark that the order of the residue integrals is not arbitrary here; the integral over 5j^i 
has to be performed first. 

The above formula is not the only possible choice for G. In fact, any other choice G' of 
the right inverse compatible with AG = ici^-, G : £i ^ £i^i can be parametrized by constants 
Aj,i,m, J, / G N, m G {1, N{1, D)}: 

G'ir" In' rhi^mix)] = G [r" In' r hi, ^x)] + 

<5f+M,,;,^r'/iz,^(x) + (5^|V2%,-r-'-^+2;^;,^(x) , (4.47) 
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where 5^ is the Kronecker delta. At each iteration step, we are free to choose one set of 

constants Aj^i^rn, ^i,Z,m- 

This freedom is partly restricted by the consistency condition. One expects that the remain- 
ing freedom corresponds to the renormalization ambiguities in the conventional framework. 
As already said, we will not prove here that our particular choice for G in eq. fl4.46p actually 
leads to a set of vertex operators that fulfill the consistency condition. 



4.2 Graphical rules for computing vertex operators 

We now want to consider in more detail the iterative scheme for calculating vertex operators 
in perturbation theory described at the beginning of this section, based on our recipe for 
inverting A that we have just described. 



4.2.1 The case D = 2: 

We will explain this first for the OPE vertex algebra in D = 2, with interaction \P{ip) = 
A ^ ^f^- What makes the construction in this theory simple is that — as we will see — the 
"counterterms" in eq. (14.411) have no effect. This is directly related to fact that the field 
theory is "super-renormalizable" in field theory parlance. As usual in perturbation theory, 
it is of no concern that the interaction polynomial is not semi-bounded, even though this 
would render impossible the non-perturbative construction of the Schwinger functions. 

In D = 2, the angular part and the A^(/, D) = 2 harmonic polynomials at order / > are 

£ = e'" , hi^^{x) = ^e^"" (4.48) 

At zeroth order, the vertex operator linear in (f is given by 

oo 

Y,{^, x) = bo In r + b+ + ^ 5^ -= [r'e^'" b+„ + r-'e^ b, J . (4.49) 

1=1 m=±l ^ 

When determining the vertex operator Yi^i{ip, x) at order i + we have to calculate Vj(<y5^, x), 
see eq. (I4.4ip . then Yi{ip^,x) and so on, and then apply the inverse G of A as in eq. (I4.39p . 
Eq. (I4.4ip contains what we called "counterterms", which are the terms with the — sign 
in front. However, it is not difficult to derive that, in D spacetime dimensions, and for a 
polynomial interaction, the scaling degree (i.e., the power of r of the most singular term 
apart from log's) of an z-th order vertex operator is 

sd {c\Yi{a, x)\b) < dim(a) + dim(6) - dim(c) + i[D - ^ deg(P)] , (4.50) 

where dim is a grading of V by the engineering dimension of the fields and is defined for a 
field of the form (Km by 

dim(a) = ^az,^(/ + ^), (4.51) 

l,m 
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i.e., each factor of ^ counts as having dimension -^j^, and each derivative as 1, which gives 
/ for the /-th order differential operator hi m{d). It follows that, in D = 2 dimension with a 
polynomial interaction, the scaling degree of all the terms Yj[ip,x — y)ip is strictly positive 
for j > 0. Hence, we may drop the counterterms in perform the limit |x — ?/| ^0 just on 
the first term. Thus, in D = 2 and i > 0, (14.411) reads 

i 

Y,{^', x) = Y, yjiv, x)Y,.^iv, x) , (4.52) 

j=0 

and an analogous equation can easily be derived for the vertex operators Yi{ipP,x). The 
only difference is that we now have p terms on the right side and a corresponding p — 1-fold 
summation. We may thus form (P' (</?), x), and hence, using the field equation (I4.38p . we 
get Yi^i{ip,x). Exphcitly, this is {i > 0) 

y,+i(^,a;) = 5^^G J2 Y,,i^,x)---Y,^icp,x). (4.53) 

For i = 0, the formula is instead [compare eq. (I3.35P 

Y,{^,x) = J2^G:Yo{^,xr: . (4.54) 

P 

Each term in the sum over p has an obvious graphical representation by a vertex with p 
incoming lines, see fig. [H 



Figure 1: Two trees representing the terms G {{G : Y^^ip , xY :) Y^^ip , x)"^) and 
G (YQ{ip,xy {G : Yq{ip,xY :)) respectively. Even though the trees are related by a reflec- 
tion, they have to be counted as different. Both make a contribution to the vertex operator 
Y2{ip,x) in the theory with interaction XP{ip) = Xip^/4\. 

Adopting such a graphical notation immediately helps one to see that if we iterate 
eq. (I4.53P starting from i = to arbitrary i, then the resulting expression will be orga- 
nized in terms of trees whose vertices have coordination number p = 1,2, . . . , deg(P), where 
P was the polynomial characterizing the interaction. We use the notation v -< w to indicate 
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that a vertex or leaf v can be reached from the vertex w following the tree downwards. By 
V ^ w, we mean {v -< w or v = w). As we have to heed the order of the creation and 
annihilation operators in the term represented by a tree, we have to count as different trees 
that are related to each other by a reflection or similar symmetry operation, see again fig. [H 
In the remainder of the subsection, we want to describe in somewhat more detail what 
the mathematical expression is for each such tree. Let T be a tree on i vertices, and let 
Yi{T,ip,x) be the contribution to Yi{(p,x) coming from that tree. Thus, 

Y,{^,x)= Yl ^^Y,{T,^,x), (4.55) 

treesronl,...,i i-i-P^' 

where is the coordination number of the vertex v. If we start from the leaves of the tree, 
then to each leaf j ending in a vertex v we have to associate a pair {lj,mj) G x {±1} 
and one of the following factors 

^ rhQi^jh'^h+^ or -^r-^^e-'^^^^^^hi.^rn- ■ (4.56) 



21 j 

Or, if = 0, we have to associate one of the factors bp or ^ ^ ^'''^bo, using the residue trick 
to generate the logarithm. The creation/annihilation operators of the leaves associated with 
the same vertex must be normal ordered, by eq. f l4.54p . It is handy to distinguish the two 
cases in eq. (14.561) by calling the lines associated with the first type of factor "incoming" , 
while calling "outgoing" the other one. For each tree, we will have to sum over all possible 
orientations of the leaves, because this will correspond to different terms contributing to 
Y.i{ip,x). 

Let us now consider a vertex v that has no further vertices attached to it downwards in the 
tree; i.e. if we follow a line downwards starting from we arrive at a leaf. At f , we have 
to multiply the factors in eq. (14.561) associated to the leaves attached to f , put this product 
into normal order and then apply the inverse Laplacian G. It is efficient to take care of the 
phase factors by introducing, for each vertex, an auxiliary integration variable < /5 < 27r, 
and to use the trivial identity 

JJ e^'".'." JJ e"''"^'^" (4.57) 

j outgoing j incoming 

= 2^/ ^ cos[/(a - /3)] JJ e''"^'^^ JJ e"'"^'^^ + (/ = term) , 

l=\ J outgoing jr incoming 

which follows immediately using that the functions form an orthonormal basis on 

[0, 27r]. The inverse of the Laplace operator has to be applied to an expression of the form 
r''cos(Z(a — /3)), where results from collecting the powers of r of the factors (I4.56P asso- 
ciated with the leaves of v. The power is thus 

E E (4.58) 

j outgoing j incoming 
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where we have assumed that none of the Z^-'s is zero. Now we introduce 5 G C \ Z and apply 
our residue trick from section ITT] to evaluate G (r'^cos(/(« — /?))). (For each incoming line 
with Ij = 0, we must replace the corresponding term with another such 6j.) The result is a 
contour integral with integrand 2 ^°'^^'g^°~^^-' r'^"^^, with u = u + 5. (This holds for / ^ 0, for 
I = we have ■^r'^^'^ .) Having introduced the additional integration parameter [3 now pays 
off, as we can carry out the sum over I using the following formula^ for G C \ Z: 



1 cos(/a;) TT cos(i^a — wk) 

1=1 ^ ' 

We can repeat this procedure for each of the remaining vertices of the tree, moving upwards 
in the tree. At each new vertex v we introduce a new integration variable (3^, and a new 
summation variable 1^ — which is summed over using the above cosine identity — as well as 
a variable defined similarly as above. More precisely, when we use the residue trick to 
perform the inverse of the Laplace operator, we must first introduce for each vertex v in 
the tree a new variable (5^, G C \ Z, and the residue in this variable then has to be taken 
in the end. If we do all this, we hence arrive at the following graphical rules for calculating 
yj(T, x), and hence Yi{ip^x): 

1. Draw the tree T with i vertices. Label the vertices by an index, f , and the lines by 
pairs of indices {vw). The vertex v has coordination number p^. The leaves also carry 
indices. 

2. With each leaf j adjacent to a vertex associate a pair {Ij, rrij) G Z4, x {±1} and one 
of the following factors 

^^m,i,(3.^+^ or ^e-'"^'^^" . (4.60) 



21 j '^'"^ ^21 j 

The first factor is chosen if the line associated with the leaf is oriented upwards, and 
the second if it is oriented downwards. For the zero modes [Ij = 0), we have to 
associate one of the factors bg or f fir^^bo, where 6j G C\Z, again depending on the 
orientation. The creation/annihilation operators of the leaves connected to the same 
vertex are to be normal ordered. 

3. With each vertex v associate a parameter 5^, G C \ Z, and a parameter f]y G [0, 27r]. 



"^This may be viewed as a degenerate case of the Dougall identity, see appendix [X] To prove the identity, 
consider the contour integral 

/ cos {z{a — tt)) dz 
J(j z{v — z) sinTTZ 

where C is a circle around the origin with radius 7r(A/ + 1/2), M € N. For M 00, the integral vanishes 
and we obtain eq. (|4.59p via the residue theorem. 
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4. With each vertex v we associate z/^ G C \ Z defined by 

E E ^^•+ E (2 + '^-)- (4.61) 

in leaves J out leaves vertices «):;<t) 

The "2" results from the inversion of the Laplacian, which at each inversion step (i.e., 
each vertex) raises the power of the radial coordinate by 2. The 5^ arises from the 
residue trick for the Laplace inversion at each vertex w below v. 

5. With the root, associate the parameter a G [0,27r], and the factor r'^-'°°^ collecting all 
the factors of r, where x = re^ "-ipha^ Th.e number i/root G C is defined as in eq. f l4.6ip . 
but with the vertex v replaced by the root, so that the sums contain contributions 
from all leaves of T. 

6. With each line {vw) (not connecting a leaf) associate a factor 

TT COs[(/3^ - (^w)^w- T^^w] 

sin(7rz/^) i/^ 
This results from the application of the cosine identity f l4.59p . 

7. Perform the sum over all Ijifaj and zero modes. Furthermore, perform the integrals 

n f n ■ 

vertices I) vertices " 

Finally, take the sum over all possible orientations of the leaves. 



The last step requires some further comment. Our first comment concerns the choice 
of the integration contours in the residue integral. They can be chosen as small circles 
around the origin with radii \5v\ chosen so that for any subset V of the vertices in the tree, 
and for any values of the phases of the 5^, f G V, we have Yliv<^v i^?; G C \ Z. We need 
this in order for our residue trick to work; the exponents Vy have to be non-integer. The 
second remark concerns the convergence of the multiple summation over the counters Ij. 
To illustrate the general point, we consider the example tree with a fixed orientation of the 
leaves given in fig. [21 

We apply the above rules and obtain the following expression for fig. [2j 



JCl JC2 "''1 ^0 ^0 , , , ^, ^„ 



h,h,h mi,m,2,m,3 



TT 



X 



X 



COS [(/5i - /32 - 7r)z/2] TT COS [{a - Pi - 7r)z/i] 



z/2 sm 7rz/2 Ui sm ttz/i 
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Figure 2: Tree contributing to Y2{(p,x), with fixed orientation of the leaves. 



with 1^2 = —h — h + '2. + S2, 1^1 = —h — I2 + h + 4: + 61 + 62, which makes a contribution to 
the vertex operator Y2{ip,x) in the theory with interaction XP{ip) = X(p^/3\. We are going 
to normal-order the product of creation and annihilation operators b^"*^ rni^t> m2^h,m3, which 
will make it easier to calculate matrix elements of the vertex operator. (Recall that the 
vertex operators were related to the OPE coefficients by the formula {c\Y{a, x)\b) = C^^(x).) 
The only terms for which normal ordering has non-trivial effects are those where either 
(/i,mi) = (/a, ma) or {12,^2) = {13,1713). These terms can be visualized by joining the leaves 
ji, or j2, is respectively by a new line, see fig. El 




Figure 3: Application of Wick's theorem to a tree. On the left side of the arrow, we have 
the tree from fig. [2l On the right side, we have the representation of the normal ordered 
expression for it. Contractions are represented by an additional line. 

We see that, in general, we will "close loops" in the tree when we apply Wick's theorem 
(see eq. flC.175p ) to the products of creation and annihilation operators. Thus, when normal 



ordering the operators in the expression Yi{T, x), we get contractions between the creation 
and annihilation operators between the individual leaves (not attached to the same vertex, 
as these already are normal ordered). 

Let j and k be leaves. Let j be outgoing and standing to the left of k, which shall be 
incoming. Their contraction gives ^i^,if.^mj,mk- We combine this with the other factors in 
eq. (14.601) and carry out the sum over rrij, getting when Ij ^ 0, where V is the 

vertex adjacent to j and w the vertex adjacent to k. When Ij = 0, we get instead 4-. We will 
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represent each such new factor by a new hne joining the respective leaves. If we apply this 
systematically, we are thus led to a wider class of graphical objects that are obtained from 
our tree T by joining an arbitrary number of leaves, but never joining two leaves from the 
same vertex. These new graphs, which we call G, are not trees any more, but also contain 
loops. The orientation of the loops is always from the "left" to the "right" , pointing from a 
leaf representing a creation operator to a leaf representing an annihilation operator. In the 
language of graphs, T is a "spanning tree" for each G, and we collect these graphs in the set 

G(T) = {graphs G \ T a. spanning tree for G}. (4.63) 

Our graphs G contain three different kinds of lines, or edges e: The edges e G T that were 
already present in the tree T, the edges e G G\T that were created by joining two leaves, and 
the leaves e that were not joined (the "external lines"). The edges in the first category carry 
momenta z/g G C that are determined by the momentum conservation rule fl4.6ip . The edges 
in the second category carry "loop" momentum l^. G Z, and the leaves e carry "external" 
quantum numbers le, rrie. The collection of loop and external momenta together is the same 
as the assignment le above. We denote by Yi{G,ip,x) the contribution to Yi{(p,x) from such 
an individual loop graph G where the sum over all possible orientations of the external lines 
is understood. In other words, 

F,(^,a;)= Yl ft^ E Y.iG,cp,x). (4.64) 

treesTonl,...,j ^^P''' G&Q{T) 

The contribution Yi{G, v?, x) from an individual graph G G Q{T), is given in turn by 

y,(G,^.x) = Y. E E f n / f'^e- 

leaf Ze,me,eleaf «e,eeG\r Vverticest; "^'^'^ ^ "^^ 
orientations 



n TI" COs[(/e + 4)(/3e - Tt)] -p-r COs(/3e/e) 

^ Sin7r(/e + 4) le + 5e h 

eST \ e ej e e e£G\T 



exp In r J J2 ~ Yl + (2 + 5, 



n b 



. in leaves e out leaves e vertices v 




(4.65) 



leaves e V e 

This formula requires several comments. First, we have used now the notation "e" for the 
edges of e, which include those present already in T, and those that close the loops, i.e. in 
G\T. To each of the latter, we have associated an integer le G N, which is summed over. 
Each of the lines e already present in T also carries an le G Z, which is determined by the 
"conservation rule" 

2= J] /e- 5^ (4.66) 

e outgoing e incoming 
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at every vertex v, where the incoming hnes e are either incoming leaves or hnes below v 
(in the sense of the relation -<), and where the outgoing lines e are either outgoing leaves 
or the line above v (in the sense of the relation >~). Hence, the numbers le,e G T are 
determined by the L loop momenta le,e G G\T and external momenta (leaves) /g via the 
momentum conservation rule at the vertices of G. That momentum conservation rule comes 
from eq. fl4.6ip . We have also introduced 6e as the sum J2v '^^ those 6y G C \ Z that 
are associated with a vertex "below" the line e in the original tree, and this accounts for 
the corresponding term in eq. 04.61 1) . Finally, also the leaves (i.e. uncontracted lines) carry 
indices le, me that are summed, and we have set j3e = Pv — if e = (vw), with = a if f 
is the root, and Pe = /3„ if e is a leaf associated to the vertex v. As always, x is related to r 
and a hj X = re'". 



w 



I 



4" 

/ I \ 



J j 



w 



1 COs[{Uw)[Pe-T^)\ 

sin(7ri/u,) 



(and impose "momentum conservation" eq. fl4.6ip ) 



]_^lj^imjljl3^ 

'21' 



]_^-ljQ-imjljl3y ^ 



COs(/3e^e) 
le 



Figure 4: The rules for the "amplitude" of a graph. The appropriate summations and 
integrals have to be understood. Slightly deviating from the earlier statement of the rules, 
we do not associate all powers of r with the root but all leaves and all vertices carry some 
power of r. 



The formula for Yi[G,(p,x) given above looks quite complicated. To see that the sums 
converge, we now introduce the following trick that, in effect, replaces the sums by integrals. 
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The basic trick is to implement the conservation rule fl4.66p by the integral 

dt 
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e=(vw):v)^w 




exp it J 2 + J2 - Yl M ' 



with one new integration variable per vertex v because momentum conservation holds at 
each vertex. Once the momentum conservation rule has been implemented (at the expense of 
the new integrals), we can now sum over all le independently, and not just the le associated 
with the loop lines e G G\T . We then do the new integrals afterwards. The advantage 
is that the sums can now be performed explicitly. We hence trade the summations over 
/e, e G G \ T for the integrals over the t^, where v runs through the vertices, as we now 
explain in more detail. 

The infinite summations that we have to deal with occur at each line e G G that is not 
an external line (leaf), and they are performed with the aid of the formula 



C0S[(/ + 5)^] U _ 1 ( ,p „ p , e . . . ^ , ^ ,.-i/3 




^ ^ ^ -A e''" 2^1(5, l;l + 6; M") + e-^*^Fi(5, 1; 1 + 5; /le"^^) | , (4.68) 



which holds for 5 G C\Z, and |/;,| < 1, and which follows straightforwardly from the definition 
of the Gauss hypergeometric series, see appendix El In the above complicated formula (14.651) 
for Yi{G, if, x), we now apply the formula at each internal line e G T, with the choice I = le, 
6 = 6e, a = Pe, and h = e'*% where = — t^, for the line e = {yw). For the loop lines 
(e G G \ T), we have the formula 

cos(/a) 



-h^ = - In + -2/icosa . (4.69) 



However, we cannot apply these formulae straightforwardly, for two reasons. First, they only 
hold for \h\ < 1, while our choice would correspond to \h\ = 1, which is on the boundary of 
the disk of convergence of the Gauss hypergeometric series. Secondly, the series in Yi{G, (p, x) 
have to go over both positive and negative values of le if e is neither a leaf nor a loop line, 
whereas eq. (14.681) only goes over non-negative values. The remedy to this difficulty is as 
follows. We first split each sum over G Z into positive and negative values of le- Then, 
for the positive values of le, we replace te by te + iO (i.e., we add a small positive imaginary 
part), and for negative values of le we replace te by te — iO. This then justifies the exchange of 
summation and integration over the t^ before we take the small imaginary part to zero, and 
the remaining integration over the ty must then be understood in the sense of distributions. 

The trick affects our formula (I4.65P for Yi{G, (p, x) as follows. Instead of the sum/integral 
in the first line of eq. (14.651) . we now have 



E E n / / ""'-d'. / • I*-™) 



leaf Ze,me,e leaf 
orientations 
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i.e., we have gotten rid of the — potentially dangerous — summation over the loop momenta, 
and replaced these by additional integrations, which are easier to control as we will see in a 
moment. Furthermore, the term in the first product (e G T) in the second line becomes 

^ in; i(te+iO)«COs[(/ + (5e)/?e] 



9,(6^, COS (3^, Q := 

sm 7ro„ V ^-^ 



sin TT^e V ^ I + 6e 



+ (4.71) 



l=—oo 



Using eq. (14.681) . we get 

g2{6e,cos(3e,te) = (4.72) 

2Fi(5e, 1; 1 + 4; -e'(+'3'=+*'=+^°)) - 2Fi(-5e, 1; 1 - 5^; -e^(-^-*^+'o)) - 1 



2Se sin 7r5e 



+ ^^r-^ f 1; 1 + -e'(-'^^+*^+^°)) - 2Fi(-5e, 1; 1 - 5e; -e^(+/^-*=+'o)) - 1 

2oe sm vroe \ 

The term in the second product (e G G \ T) is as in eq. (I4.69p . The remaining parts of 
the formula (I4.65P are unchanged. The total effect of these manipulations is reflected in 
eq. (15.114p . setting D = 2 and e„ = there. 

A possible divergence can now only come from the d(3ydtv integrations, and the danger can 
only come from configurations near f3e = ±te + 7r, where the argument of the hypergeometric 
function tends to 1. It is at this stage that having integrals instead of sums pays off, because 
we can now use the well-known expansion formula for the Gauss hypergeometric function 
near x = 1: 

2F1 (5, 1; 1 + 5; x) = 26 [7e - ^(5)] - 1 - 5 ln(l - x) + o(l - a;) . (4.73) 
When we apply this to the 2-^1-factors in the terms above, we find that 

g2{Se,COSpe,te) ^ . . . ^ ^ lu - - 7r| + In \pe + te - 7l\ , (4.74) 

near /3e = ±te + tt, where we have also used standard identities such as ln(x + iO) = In |x| + 
i7r9{x). The contour integrals over the parameters 6^ are harmless as long as we choose 
the contours such that dist(5e,^) > const. These estimates suffice to show that any matrix 
element {c\Yi{ip,x)\b) is convergent. To see this, we recall that Yi{ip,x) was the sum over 
all trees T on i elements of the quantities Yi{T,ip,x), which in turn was the sum over all 
G G QiT) of the quantities Yi{G,ip,x), for which we gave a formula above. Now, when 
forming a matrix element {c\Yi{G, (p, x)\b), we get a sum of terms of the form 

n K,rr..^^-- b). (4.75) 
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We write (c| = {c\)~^/^{0\ U (b;,^)''-'" and \b) = H (b^^)'"''" |0). Furthermore, for a 

fixed assignment of orientations and indices le, rrie to each of the (uncontracted) leaves e G G, 
we introduce the multiindices a~^,a~ by 

(^tm = # {e e G : e incoming, {le, rrie) = {I, m)} 

%m = ih{(^^G : e outgoing, (/e,"^e) = . (4.76) 

In order for the matrix element (I4.75P not to vanish, all creation operators have to be con- 
tracted with annihilation operators of the same indices. Remember that we are considering 
the contribution of a graph G with a fixed number of leaves. In the above notation, that 
means \a^\ + \a~\ is fixed. We also assume that the multiindices 6, c are fixed. This means 
there is maximally one choice for a+,a~ so that the matrix element (14.751) does not vanish. 
In this case the latter is equal to 

623= n -7s^ n -7^- (4-77) 



/2/e M 

e incoming e outgoing 



e 



The sum over the orientations and indices of the leaves in eq. (14.651) thus reduces to the 
sum over all graphs with la"*"! incoming leaves carrying indices from a+ and \a~\ outgoing 
leaves carrying indices from a~ . This sum is obviously finite. 

Thus, the only source of divergent behavior are the integrals over the f3v,tv By our 
estimates above, these are bounded from above by (a constant times) the integral of the 
form ^ ^ 

n (f^d(3, rdAl[\n\te±Pe-7i\. (4.78) 

verticest;inG ^"^^ "^^ ^ eSG 

Such integrals are standard "loop integrals" in quantum field theory (in position space), 
and they are well-known to be finite. Hence we conclude that our formula for Yi{ip, x) gives 
a finite answer, despite the infinite sums. 



4.2.2 The case D>2 

Let us now consider the field theory with interaction XP{ip) = X ^^p^, in dimension D > 2. 
The main difference to = 2 is that eq. (I4.52p . and hence eq. (14.531) . no longer holds, as 
the "counterterms" in eq. (14.411) now do not have a vanishing limit as |x — ?/| 0. This 
is closely related to the fact that, for D = 2, we were dealing with the special case of a 
super-renormalizable theory whose divergences are of a particularly benevolent nature. In 
dimensions D > 2, we must work with a version of eq. (14.411) that also incorporates the 
counterterms, and this will lead to correspondingly more complicated graphical rules. 
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Our aim is to express the vertex operator Yi{(f, function of the aheady known 

0-th order vertex operators, and then find a graphical representation for this expression. 
We hst the equations that we need to decompose Yi{ip,x) into 0-th order operators. Eq. 
fICTD reads 

Yi+,{^,y) = ^GY,{^\y). (4.79) 
Moreover, we need the equations fl4.4ip and 

i i ^ 

(4.80) 

where the dots stand again for terms vanishing in the hmit x — y ^ 

We do not perform the hmit x — y ^ for the moment, so each time we use either of the 
eqs. f l4.4ip or fl4.80p . we have to introduce a new variable from MP . We can choose this new 
variable (x above) to be (1 + e) times the old variable {y above, i.e. x = {1 + ^)y), where a 
new regulator e^, > has to be introduced each time we apply either of eqs. (14.411) or (I4.80p . 
The result of repeatedly applying eqs. (14.791) . (14.411) and (I4.80p is a sum of products of nested 

0- th order vertex operators whose arguments from depend on the initial variable from 

and the e^j's. We now focus on a special partial sum, the sum of "tree-like" summands. 
We call a summand tree-like if, when tracing back its path through the iteration process, 
at each iteration step it does neither belong to the counterterms nor to the smooth terms 
represented by dots. Another way to put this is to say that we are only interested in the 
terms that we would have obtained by dropping the counterterms and dots in eqs. (14.4 ip 
and fICTD . 

A formal definition of the tree-like summands can be given as follows. Let T be a tree 
such that all vertices have coordination number four. With each vertex or leaf v of the tree 
associate a regulator e^. If we have three such trees Ti,T2,T3 with number vertices ii,i2,i3, 
then we can join these at their root to form a new tree, T = U^T^^ with i vertices. The three 
vertices now attached to the new root (corresponding to the 3 subtrees Ti, T2, T3) carry new 
regulators denoted e.g. ei, 62, £3. The recursive definition of the tree-like contributions to the 

1- th order vertex operators is then 

1 ^ 

y,(T,<^,x) = - J2 lly^AT^,^,Xw), (4.81) 

where x^ = (1 + eu,)x, and where the dependency of the quantities Yi(T,ip,x) etc. on 
the regulators {e^} has been suppressed to lighten the notation. The vertex operators 
Yi{ip,x) are given by a sum of these tree-like terms [see eq. (I4.55p ]. plus counterterms. 
Moreover, it will turn out that these tree-like summands can be thought of as building 
blocks for the (more complicated) counterterms. We will not discuss these here, but 
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outline their construction in the next section. The regulators {e^} cannot be sent to zero 
before taking the sum of all contributions to the vertex operator, including the counterterms. 

We would now like to find a closed form expression for the tree-like terms, applying 
the same kind of reasoning as for D = 2. We have to take into account some differences. 
The first difference is that the trigonometric functions in the expression for the free vertex 
operator are now replaced by the harmonic polynomials, see eq. fl3.28p . The analogue of 
relation fl4.57p is 

n ''h^mM n (4.82) 
j outgoing j incoming 

°°2/+n — 2/" — 

/=0 D S j outgoing j incoming 

using this time the orthogonality of the harmonic polynomials, as well as the formula for 
the Legendre polynomials P{z, I, D) in D dimensions given in appendix |Al We proceed as 
in the case of D = 2 dimensions. Let us now consider a vertex v that has no further vertices 
attached to it downwards in the tree. Let the leaves attached to v be labeled by j. We collect 
corresponding factors of r, and the harmonic polynomials. The harmonic polynomials are 
multiplied by the formula just given, while the factors of r work out as r'^", where z/„ is now 
given by 

E E (^. + ^-2)- (4.83) 

j incoming j outgoing 

Thus, in total we get a factor of {21 + D ~ 2)r^^ P(x ■ /, D). We have to apply the inverse 
of the Laplacian to this expression using the residue trick, introducing 5y and z/^ = z/^, + 5^. 
When we do this, we get a contour integral with integrand 

{2l + D-2)ni-y,l,D)r'^' 

v,{v, + D-2)-l{l + D-2) ■ ^ ' ' 

The sum over I is now readily performecjf] using the generalized Dougall's identity^ see ap- 
pendix [HI which states that for any z/ e C \ Z and —1 < z < +1 and > 3, we have the 
identity 

t (2i + 0-2)P(.;i.g) 

j-^^v{v + D-2)-l{l + D-2) sinvrz/ ^ ' ' ^ ^ ' 

From the vertices connected to the leaves, we then work our way upwards, repeating for 
each new vertex the same procedure. We will then end up with a similar set of graphical 
rules. The main difference to Z) = 2 is that we have to use Legendre functions instead of the 



^This key observation is due to [H 
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trigonometric functions, and that we must take into account the regulators e^. As one can 
see, this means that we have to introduce an extra factor 

vertices v>zj 

(X -lj-D+2 
(1 + 6,) n (1 + ^-)) (4-86) 
vertices v>zj / 

for each leaf j, depending on whether it is incoming or outgoing, to get the correct rules for 
D>2. 

To summarize, we have the following graphical rules for calculating Yi{T, ip, x), the regularized 
contribution of a tree T to the vertex operator Yi[ip,x) for the theory with interaction 
XP{(p) = X(p^/4\ in D > 2 (a general interaction polynomial is completely analogous): 

1') Draw all trees with i vertices with incidence number 4. Label the vertices by an index 
V and the lines by pairs of indices (vw). The leaves also carry indices. 

2') With each vertex v associate a parameter 5„ e C \ Z, a parameter y„ e S^~^ and a 
regulator e„ > 0. 

3') With each leaf j adjacent to a vertex v, associate a pair {Ij, m) e N x {1, N{D, I)} 
and one of the following factors 



hlj,mj{yv)h+r. 



The first factor is chosen if the line associated with the leaf is oriented upwards, and 
the second if it is oriented downwards. The creation/annihilation operators of the 
leaves connected to the same vertex are to be normal ordered. Also, write down the 
regularizing factor 

((l + ejOrivertices^tjll + e^))^' 



4') With each vertex v we associate Uy e C\Z defined by 



or ((l + 6,)nvertices.^,(l + ^.))"'^"''^' (4-88) 



^v= E (^i + ^-2)+ E (2 + 5J. (4.89) 

in leaves j'^D out leaves j-<i) vertices ui:;<t) 

The "2" results from the inversion of the Laplacian, which at each inversion step (i.e., 
each vertex) raises the power of the radial coordinate by 2. The 5^ arises from the 
residue trick for the Laplace inversion at each vertex w below v. 
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5') With the root, associate the parameter x G S^^^^ and the factor r'^™"*, where x = rx. 
The number i/j-oot € C is defined as in eq. (14.6 ip . but with the vertex v replaced by the 
root, so that the sums contain contributions from all leaves of T. 

6') With each line (vw) connecting vertices v,w associate a factor 

71 

Pi-Vv ■ Vw] i^w, D). 

sm 7111 w 

This results from the application of the Dougall formula. 
7') Perform the sum over all lj,fnj. Furthermore, perform the integrals 

n I my.) UirJ/f- 

vertices v vertices v " 

Finally, take the sum over all possible orientations of the leaves. 

If we proceed as in the case D = 2 and perform the sum of all tree-like terms as in 
eqs. (14.641) . (I4.55P we do not get the complete vertex operator Yi{{p,x), so in this sense our 
rules are not complete. To get the complete vertex operator, we should also incorporate 
the counterterms in D > 2 (see next section). The sum over graphs (I4.64p depends on the 
regulators e^, and is divergent for 0. The expectation is that these divergences are 

canceled by counterterms. We will briefiy expose this idea in the next subsection, without 
explaining it in full. 

Let us write down explicitly the contributions Yi{G,(f,x) from loop graphs G as in 
eq. (I4.64p . In D dimensions, the formula is 

Yi{G,Lp,x) = 

X 



X 



X 



X 



leaf (e,me,eleaf le,e£G\T \verticost) 
orientations 



eGT ^ eeG\T 

expflnrJ " Yl i^e + D - 2) + ^ (2 + 5,) i j 

\ ^. in leaves e out leaves e vertices v ) / 

n fnd+'jV' n iw 



le / \ -le-D+2 



inleavese \v>-e / out leaves e \v)~e 



leaves e 
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In the last line, we have set jje = Vv for an uncontracted leaf e associated to some vertex v. 
The Of. are defined as follows: If f is a vertex, let 

= In J] (1 + e^) , (4.91) 

and if e = (vw) E G\T, let 9e = 9^ — 6w The regulators are always chosen so that 9e > 0, 
but evidently, ^ as the regulators are removed. 



w 



I 



sin m/^ 



4" 

/ I \ 



V 



— > r^"'"'''' (and impose "momentum conservation" eq. fl4.6ip ) 



I 



V 




Figure 5: The rules for the "amplitude" of a graph in dimension D > The appropriate 
summations and integrals have to be understood. As in fig. HI we associate the powers of r 
with the leaves and vertices instead of with the root. 

It is possible to make the source of divergences more transparent by introducing the 
same trick as in D = 2 to replace the summations over the angular momentum numbers 
le by integrals. At each vertex, we have the momentum conservation rule eq. (14.661) . As 
in D = 2, we can realize this rule by introducing an integral similar to (14.671) at each such 
vertex. We can then perform the summations over the 1^ as we have done in eq. (14. 72 p . 
The overall result of these manipulations is the following. Instead of the sum/integral in 
eq. (14.901) . we now have 



leaf Ze,me,eleaf \verticesti 
orientations 
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i.e., we have gotten rid of the summation over the loop momenta, and replaced these by 
additional integrations over the parameters The term in the first product (e G T) in the 
second line of eq. (14.901) becomes, with /3e := arccos(y„ ■ yw) for the line e = (vw) connecting 
the vertices v and w: 

oo 

gD{Se,COS(3,,Q := -iL^(V(-l)V(*^+'°)'P(cOS/3e,/ + 5e,/^) 

^ 1=0 

+ J2 (-l)'e^(*'=-'°)'P(cos/5e, I + 5e, D)) . (4.93) 



l=—oo 



Putting together eq. (14.721) and eq. ( IB. 1461) . we get for even D 



g-it,(D/2-l) / g \-V 

r(L'/2) smnde \2smpeOPe/ 
X ^e+^'^^^^ /(4, 1 + ei(+/5^+*=+io)) + e-^'^^^^ f{-6e, 1 + e'(-'3'=-*'=+'°)) 

+ e-'^=^^ /(4, 1 + e'(-^^+*=+'°)) + e-^^^^-^ /(-4, 1 + e^(+/5-*'=+io)) - cos(4/?e) 

Here, / is the transcendental function given by 

f{5,x) = y"^[ij{n + l)-ij{S + n)-\nx]x'', (4.95) 
n! 

n=0 



also see eq. flB.144p . For odd D, there is a similar formula obtained with the aid of theorem [2] 
of appendix [Bl 

The term in the second product (e G G \ T) in eq. (14.901) is affected by the regulators. 
For these edges, we must replace the term in the second product over (e G G \ T) by the 
expression 

^e"^(*=+'^^)p(cos/?e,/e,I^) = [(e^'^^ -e^(*'=+^^=))(e-'^= -e'(*=+^''^))]"^''''^^' . (4.96) 

There is also an additional factor of Y[v{vw)=e(^G\T ^^^"^^'"^^ arising from the fact that the 
momentum conservation rule in D dimensions is slightly different. The remaining parts of 
the formula (I4.90p are unchanged. The total effect of these manipulations to eq. (14.901) is 
reflected in eq. (I5.114p below. 

The divergence that appears when we set e^, = in our graphical rules for the tree-like 
terms Yi{T, (p, x) now manifests itself as a pole of the transcendental function f{6, x) a.t x = 0, 
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which corresponds to (3^ = =t^e — vr in the above formula fl4.94p . More precisely, using the 
expansion of the function / above, we can infer that the divergent behavior oi go are sums 
of terms of the form ~ (te ± /^e — vr ± iO)^~^. It can be seen from this that the most divergent 
part of the integrals that we have to consider is now given by an expression of the form 

Ig= n il^, dnAt^\{\te±Pe-A''''' W \te + iee±Pe?-'' . (4.97) 
vertices in G S'^-^xS^ / ^^rp eeG\T 

In D > 3, integrals of this type are no longer convergent for e^, = (i.e. 6e = 0), but the 
divergences are very similar in nature to the divergences found in ordinary Feynman integrals 
in x-space. (But note that the "x" space has become S^^^ x S^, and that the iO prescriptions 
are different). In the usual approach to perturbation theory via Feynman integrals, these 
divergences have to be "renormalized" by hand. But in our approach the divergences are 
expected to cancel automatically when we include all contributions to Yi{ip, x) including the 
tree like terms and the counter terms before taking the regulators to zero. We now explain 
this in a little more detail. 



4.3 Renormalization 

In conventional perturbative quantum field theory, renormalization is necessary to make the 
terms in the perturbation series well-defined. In our approach the situation is somewhat 
different, because the consistency condition tells us in principle right from the start how to 
obtain well-defined perturbations of arbitrary order. We think this is a remarkable feature 
of the present framework. Nevertheless, we have already borrowed some vocabulary from 
renormalization such as "counterterms" , and the reason for this is that we also need to 
perform various limits in our approach which are quite reminiscent of certain operations in 
conventional renormalization theory. In fact, the inclusion of the counterterms into the rules 
may be thought of as the "renormalization" of the tree-like contributions that we have treated 
in the last subsection and that diverge when the regulators are sent to 0. The counterterms 
cure these divergences. 

To state the "complete" set of rules is in principle straightforward, but generates a rather 
heavy notation. To keep the discussion reasonably transparent, we therefore only sketch the 
basic procedure, and leave a more detailed discussion to a future paper, where we will also 
outline the relation to Hopf-algebras similar to those in [7]. 

Let us go back to the start of the recursion procedure for evaluating Yi{{p,x). By 
eq. drZS]), 

Y,{^,x) = ^GY,.,{^^x). (4.98) 
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By eq. (giHOD, 

i-l 
j=0 

- J2Yi^,^,{Yj{ip,ex)ip\x)---^Yi.-i{ip,x) + ... (4.99) 

where dots stand for terms vanishing for e — > 0. The counterterms (the terms with the 
— sign in front) can be rewritten as 

^ Yi_i^j(Yj{(p, ea;)v3^ x) + - — r^Yi_i{ip, x) 
i=i 

i-i ^ 
= Yl Yl (c,F,(^,ex)¥p2)y,_i_,(c,x) + -^^F,_i((^,x) + ... (4.100) 

j=l dim(c)<3{D"2)/2 ' ' 

where again, dots stand for terms vanishing for e — 0. None of the terms in eq. f l4.100p will 
give rise to tree-like summands in the final formula for Yi{(p, x). However, we can now apply 
the "incomplete" rules from the last subsection to the operators Yj{(p,ex), Fj_i_j(c, x) and 
Yi_i{ip,x) appearing in eq. (14.1001) . and thus obtain more terms that contribute to Yi{(p,x). 
Graphically, we represent this by a "decorated tree" where the trees contributing to the 
matrix element {c\Yj{(p,ex)\(p'^) are the decorations of a special vertex or blob, to which in 
turn the roots of the trees contributing to Fj_i_j(c, x) are attachecj^, see fig. El 

Of course we still do not have identified all the terms that make a contribution to Yi{ip, x). 
We would have to take into account the counterterms at each recursion step in the same 
manner as above. But we see that if we do so, we get all summands in the final regularized 
formula for Yi[ip,x) represented by (multiply) decorated trees. To summarize, we have laid 
out the following idea for calculating a vertex operator Yi{ip,x): 

• Draw a certain set of decorated trees with i vertices. 

• Apply a set of rules similar to those from subsec. l4.2.2l to each of these trees, translating 
them into an "amplitude" depending on a number of regulators e^. 

• Take the sum of all those amplitudes and take the limits ^ in the appropriate 
order. 

The development of rules for drawing decorated forests mentioned in the first bullet point 
above and the rules for amplitudes of decorated trees from the second bullet point is 

^The term \€x\'^^^Yi^i{ip, x) would be represented in the same way by a blob decorated with the matrix 
element |ea;p~^ — {(f\Yo{€x)\ip'^) , to which a tree contributing to Yi-i{ip,x) is attached. 
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Figure 6: A decorated tree making a contribution to the vertex operator Y4^{ip,x). More 
precisely, this tree makes a contribution to G {{ip\Y2{ip,ex)\ip'^)Yi{(p,x)); the tree pointed at 
by the arrow is a summand to Y2{(p,ex) and the tree below the "blob" is a summand to 
Yi{ip,x). 



straightforward but cumbersome. They naturally lead to the appearance of a Hopf algebra 
structure similar to that found in [7] for ordinary Feynman diagrams. We will not state 
these rules here, leaving this to a future paper. 

As an example, we consider the operator Y2{<f,x) for the theory with interaction 
X(f^/3\ in D = 4. There are two trees with two 3-valent vertices. We take the sum over all 
graphs spanned by these trees and include the one decorated tree that makes a contribution 
to the sum, see fig. [71 




Figure 7: All trees contributing to the operator Y2{ip,x) for D = 4, P{f) = (p^/3\. In the 
decorated tree, the tree pointed at by the arrow is contracted with the vectors {(p\, \(p). 
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The contribution of the decorated tree is 



52(2 + 52) 4/ + 2/((5i + ^2) + 8 + 6(5i + ^2) + (5i + ^2) 



(4.101) 



Here, we have aheady performed the (trivial) integrals over the angular variables and we 
have not used the trick fl4.85p because the sums over the "momenta" of internal lines are 
trivial as well. In the limit e ^ 0, the expression (14.1011) is logarithmically divergent, as 
can be seen from evaluating the contour integral in ^2. This has to be compared with the 
divergences in the loop graphs. Consider the contribution from the second tree in the first 
line of fig. [3 



Ci "1 JC2 "2 Js'i 



2 



X P{y, ■ ^2, k, 4)-^ P(yi -^,1-11 + 62, 4)(1 + ey-'^+''' 
sm 7ro2 

^ J , . , P(yi-x,/ + 2 + 5i + 52,4) 

sm7r(di + ^2) 

X r'+2+^^+^^/i,,„(y2)cu(4,/)-^/2b+^ (4.102) 

We have only kept track of one regulator e; the other regulators that one would obtain 
by following the graphical rules by the letter give rise to trivial limits e„ — > 0. 

The contributions from the other loop graphs are similar. In the loop graphs, the 
divergences arise from the loop summations which is the sum over li in eq. (14.1021) . These 
divergences should be canceled by the counterterm (14.1011) . In the previous section, we 
have applied manipulations which turned the sums into integrals, and the divergences then 
manifest themselves at the level of the integrals Iq, as we discussed. It is still not an easy 
matter to see that the counterterms precisely cancel the divergence, already in this example. 
See [13] for a detailed discussion of examples of that nature. It seems that a different 
method is required to make such cancellations manifest in a systematic way. 

Counterterms also occur in the construction of composite operators. Here, what have been 
"tree-like summands" before become "forest-like summands". For an operator Yi{a,x) these 
would have a graphical representation by forests consisting of |a| trees with i vertices overall. 
The rules for the contribution to Yi{a,x) of such a forest are more or less the same as the 
rules from section 14.2.21 applied to every tree in the forest. Again, the notation necessary to 
state the rules in full, including decorated forests for the counterterms, is quite heavy, and 
we will not spell this out here. 



37 



5 Vertex algebra and special functions 

In this section, we will give yet another representation of the vertex operator Y{ip,x) in the 
interacting field theory characterized by the interaction A P{<f) = A ^ ^f^- believe that 
this representation is interesting not only because it might be convenient for calculations, 
but also because it hints at a deep relation between our vertex algebras and certain special 
functions of hypergeometric type. 

Our starting point is the representation of Y{ip,x) = A*yi((/9, x) in the form of 
eqs. (14.651) and (14.901) . We further decompose the i-th order contribution Yi{ip,x) to the 
vertex operator into contributions from counterterms (when D > 2), and from graphs, G, 
as outlined above, and these are called again Yi{G,ip,x). The main point of this section 
is to provide an alternative representation for this quantity. Since each Yi{G,(p,x) is an 
endomorphism in End(K), it is sufficient to present a formula for the matrix elements in a 
basis of V. 

To present our result in the most economical form, it is convenient to choose a particular 
(non-orthogonal) basis of V which is defined as follows. First, for any p G / G N, we define 
hiipy = Lu{l, DY^I^ hmip) Km- For p = (pi, . . . , p J G M"^ and 1 = (h, ...,/„) G N", 
we then define 

n 

\p,r):=l[hiM\0). (5.103) 

i=l 

We remember that any element in V corresponds to a composite field, i.e., a formal product 
of ip and its derivatives. For the vectors just defined, this is 

n 

\p,r) = l[P{p,-d,h,D)ip. (5.104) 

1=1 

It is evident from this expression that the vectors \p, l) form an (overcomplete) basis of 
V. Our aim in this section is to provide an alternative expression for the matrix elements 
{p^,l^\Yi{G,(p,x)\p+,l+). 

Our starting point is the formula (I4.90p for the Yi{G,(p,x). The basic idea is to carry 
out the angular integrals J df2(y^) first, or rather, to turn these integrals into summations. 
This is done by first expanding the Legendre functions using the formula (valid for z/ ^ Z) 

which we prove in the appendix [Bl For z/ = / G Z, we have instead the well-known formula 

2-(D-3)/^m r(l-j + D/2-l) 
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If we perform these expansions for all the Legendre functions appearing in eq. fl4.90p . then 
we end up with a multiple sum, whose terms contain powers {y^ ■ where e = (vw) is 

an edge between v,w, and where each rig is the summation counter from the power series 
expansion of the Legendre functions associated with the edge e. To perform the angular 
integrals, we now further expand each such power using the multinomial formula. 



Here, and in the following, runs from 1 to D. After the combined expansions, each term 
in eq. (14.901) will now consist of a prefactor times yv,^^■, raised to some power a„^^. The power 
is 

a^,/. = E ^'^'A" (5.108) 

e adjacent v 

where the sum is over all edges e going from the vertex v. Thus, the integrals we have to 
consider are of the type (aj G N) 

2^tI^-^— if 6ven , 
dfi(x)xf ■■■x^^ = <( n^^^) (5.109) 



SD-l 



otherwise . 



This formula can be viewed as a multi-dimensional generalization of the standard formula 
for the Euler Beta-function {D = 2) and can be proved e.g. by induction in D, expressing 
dQ in Z)- dimensional polar coordinates. If we combine all the steps we have described so far, 
then we end up with the following expression (up to a numerical prefactor) for the vertex 
operator: 



{p^,L\Y,{G,if,x)\p+,l+)= (5.110) 

1 r d6: 

27ri Jc, 



^2. 



n 



E E E n^/„^ 

assignments l^eN: e(^G\T ke&°:eeT \veT 
i+,i_—> leaves 

r(-/e/2 - 6J2 + \k,\/2)T{lj2 + 6, + D/2-l + \k,\/2) 



kj 

l-r n^r((Eeon.fee,^ + l)/2) yr T {k - J, + D /2 - 1) 

yr r(/+e - Je + D/2 - 1) yr TjU ' Je + D/2 - 1) 
M m^e-2j^l 11 Je!(/-e-2,e)! 

X ffc0^Eein« + e-E.outa-e-^ + 2) + E„eT^- (_2)^-l'''=l JJp^'^g JJp^'^g H^"'' 

ein eout e£G\T 
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Eq. (IS.llip requires several comments. First, we have the regularizing factors from eq. fl4.90p . 
This means that eq. (15.1111) rigourously makes sense in arbitrary but for D > "i the 
expression is divergent if we let — > 0. We have left out the regularizing factors associated 
with uncontracted leaves because they remain regular in the limit e„ — >■ 0. As above, T is 
a tree on 1, . . . , z, and G is a graph from Q{T), the set of all loop graphs for which T is a 
spanning tree, and which result from T by joining leaves together and hence forming loops. 
The edges are e, and the vertices v. The leaves on T that were joined to form loops are 
hence e & G \ T, and with each of them, we have a summation counter /g G N. The l^s 
associated with the uncontracted leaves are assigned to the entries in l±. Here, the notation 
"e in" means that e runs through the incoming leaves, while "e out" that it runs through 
the outgoing leaves. The le associated with the incoming leaves is set equal to an Z+e in 
the vector /_|_, while the Ig associated with the outgoing leaves is set equal to an /_e in the 
vector Z_. By the first summation sign in the first line of eq. f lS.llip . we mean the sum over 
all possible assignments of the entries of /+,/_ to the uncontracted leaves of G. A similar 
notation is used above for the vectors p+,p-. 

The le associated with edges that are neither leaves, nor contracted leaves are referred to 
as "e G T", and if e = (vw) they are also called in the last line. The value of le on such 
an edge is determined by the momentum conservation rule f l4.66p . Furthermore, with each 
edge, we have additional summation counters ke G N^. The sum over ke when e & G \ T 
is over those integer vectors for which 2je := le — \ke\ is even. For such integer vectors 
k = {kfj) G N-^, and vectors p = (p^) G M.^ we use the following multi- index notation: 

\k\ = J2K, P' = l[pt'^ k\ = l[k,\. (5.111) 

H n 

We write "e on v" to mean that the sum/product is running over those edges e going out 
from the vertex v. The integer vector ko G is the counter associated with the root e = 0. 

Formula (15.1111) is the desired alternative representation of the contribution to the vertex 
operator from an individual graph G. In D = 2, the full vertex operator is given by the sum 
of the graph contributions (15.1111) with = 0, as in eq. (I4.64p . For D > 3, we have to keep 
the regulators finite. The full vertex operator is the sum of the contributions from the 
graphs (I5.11ip and the counterterms whose construction was outlined in sec. 14.31 The limit 
as has to be performed in the end. 

The residue integrals J 6.6^/6^ can be performed straightforwardly using the well-known 
Laurent expansion of the Gamma-function around integer values, which can be inferred from 
the standard formula 

r(l + 6) = ^e^(^-^-)exp |f^(-5)"(Cn - l)/r^| , (5.112) 

where (n are the values of the Riemann Zeta-functioiill. Thus, we see that we get a repre- 
^We expect to see a connection to the "Z-sums" described e.g. in [37] and their algebra when performing 
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sentation involving only (multiple) sums. 

It is worth contrasting the above formula flS.llip with the integral representations for 
the contribution Yi{G,ip,x) that were derived earlier in sec. HI For comparison, we here 
repeat these formula, or more precisely, the matrix elements of them. In doing so, we 
take the opportunity to replace the integrals Ylv Is^-i ^^{Vv) over the spheres by integrals 
over the corresponding 5'0(L')-invariants ■ y^, using a general formula due to [TU]. More 
precisely, let us write {p+,'P-} = {pi \ i = I, L} for the "momenta" associated with 
the states \'p±,l±). We order the vertices so that the vertex f = is the root, so that 
the vertices v = 1, . . . , L correspond to the leaves, and so that the internal vertices are 
v = L + l,...,L + i + l. Then, for each pair (vw) of vertices we introduce a real integration 
variable Zy^^, and if (vw) = e represents an edge e of the graph G, we also write z^^^ = Ze- 
From these quantities, we also define the following symmetric D x D matrices: 



ZD-2,0 



Zk 



\ 



\ Zo^k 



ZD-2,D-2 ^fc,D-2 
ZD-2,k Zk,k J 



(5.113) 



The formulae derived in sec. Hlcan then be seen to take the form m D > 2 (up to numerical 
factors) : 



X 



X 



X 



X 



X 



(p_,L|ri(G,v2,a;)|p+,/+) = r 



2i+E.in' + e-Eeout(«-e + ^-2) 



n 



assignments \L+l<.v<.L+i+l 
Ij^ .(—-^leaves 

r2-K 

e'^^-dt. 



1 

27ri 



n 



l<v<L+i+l 







n 



l<v<L 



ZO,: 



■^D,L+^+l L+i<v<L+i+i \/det(Zz), 

n 



0<v<D-2 
v<w<L+i+l 



l<v<D~2 
v<w<L 



Y[gD{Se,Ze,te) JJ [(& 

eST eeG\T 



1 arccos Ze 



gi(ie+i6'e)'^^g-i arccos Ze _ g 



l[F{z,,l+,,D)l[F{z,,U,D) n 

ein eout l<f<i'+ 



i(ie+i^?e)^l-(^-2)/2 

(5.114) 



Here, qd are the distributions given explicitly in eq. (14.941) . and the (5's in the fourth line are 
Dirac 5-distributions. In D = 2, the terms in the last product are replaced by log's. 



the residue in eq. (jS.llip . 
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Furthermore, ^D,n is the manifold of dimension (r — l)(n — r/2) of all real, positive n x n 
matrices M = (Mij) of rank r = min{D,n) such that Ma = 1. Thus, we see that while 
formula flS.llip is given entirely in terms of infinite sums, the above representation fl5.114p is 
entirely in terms of convergent integrals. The integral over ^0,1+1+1 will become divergent 
for D > 3 in the limit when the regulators are removed, the divergences coming from the 
boundary d^n^L^i^i. As was the case for the previous formula, for D > 3, we would have 
to calculate the vertex operator {p^,l-\Yi{(f,x)\p+,l+) as the sum of all the above terms 

/+) plus the counterterms, before removing the regulators. This was 

outlined in sec. 14.31 

It is fair to ask what is the value of having the alternative representations (15.1111) 
and (15.1141) . It is not clear that either representation has much of an advantage computa- 
tionally, as there is essentially an equivalent number of summations as there are integrations 
in both formulae (15.1111) and (I5.114p . However, the alternative representation (I5.11ip brings 
out a striking feature that was far from obvious when we started the construction of the 
vertex operators, namely that it can be represented in terms of (multiple) infinite series of 
a very special form, with each term being a monomial times a ratio of Gamma-functions. 
Because of this feature, the above series can be viewed as a generalization of the Gauss hy- 
pergeometric series, and the particular form of the series is governed by the graph G under 
consideration. The vertex operators Yi{a, x) may also be defined for general a & V, and their 
matrix elements have a similar representation. Furthermore, the vertex operators satisfy the 
consistency relation (12. 9p . We expect that this relation will give highly non-trivial relations 
between the above functions of hypergeometric type. Those relations make this class of 
multi-variate functions special. By analyzing the relations that are obtained in more detail, 
we expect that one can uncover interesting relations between our vertex operator algebras 
and the theory of special functions. We think that this is a fruitful direction for further 
research. 

It remains to be understood better how to incorporate the renormalization procedure into 
this approach. As we have emphasized before, the above formula only yields a convergent 
result for D = 2, while it is only an incomplete representation of the vertex operator for 
D > 2. It appears that a formula of the above nature including renormalization can also be 
given, but this involves the use of a Hopf-algebra structure similar to that of [7]. We will 
pursue this in another paper as well. 

6 Relation to vertex operator algebras in CFT's 

Vertex operator algebras have been discussed previously in the literature (see e.g. [3l [Qj [T9]) 
in the context of conformally invariant QFT's (CFT's), so we explain the difference between 
our framework/motivation and those approaches. 

First of all, we stress again our framework is in essence a "repackaging" of the information 
contained in the Wilson operator product expansion (OPE), which is a standard tool in 
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QFT. The relation is, as we explained above, {c\Y{a,x)\b) = C^f^{x). Our new realization 
is that, when repackaged in this way, the OPE satisfies properties that can be encoded in 
some kind of vertex algebra. As we have discussed, this repackaging allows both for a new 
conceptual viewpoint of perturbation theory via Hochschild cohomology and a new way to 
do calculations in perturbative QFT. 

By contrast to the vertex operator algebras that have previously been considered in the 
context of conformally invariant theories in D = 2 dimensions, our approach is intended to 
work not just for 2-dimensional CFT's, but for any QFT (in any dimension) whose Schwinger 
functions have an OPE. Despite this key difference, there are some evident parallels: As in 
our approach, one considers endomorphisms Y{a,x) G End(V^), where a eV. Also, similar 
to our approach, the vectors a are interpreted both as states and as fields ("state-field 
correspondence"), and the Y{a,x) satisfy certain properties that are similar to ours. More 
precisely, in the usual approaches, is identified with C, and the vertex operators are 
formal distribution^, Y{a, . ) G End(l^) C[[x, x~^]], where C[[x, x~^]] is the ring of formal 
sums of the form X^fcez ^-fca^^ with complex coefficients. The analogue of the transformation 
formula under 5*0(2) is played by a formula expressing the conformal invariance of the 
theory, i.e. either under PSL2{M.) or even the Virasoro algebra. A notion of vertex algebra 
in D dimensions that is in a similar spirit and is applicable to globally conformally invariant 
theories has also been introduced, see [2^. The appearance of the ring of formal Laurent 
series in the CFT context is connected in an essential way to the fact that, in conformally 
invariant theories, the OPE-coefficients only have singularities of a very special form due to 
conformal invariance. 

Our consistency condition, as stated above in eq. (11.51) . does not make sense as it stands 
in the standard CFT- vertex algebra context P, [19], because the left side of the condition 
would be an element in the ring End(K) ® C[[a;, x~^,y, y~^]], whereas the right side would be 
an element in the ring End{V) (g) C[[{x -y),{x- y)~%[y, y'^]], and there is no natural way 
to identify these rings without some notion of convergence. Such a notion of convergence is 
available (and used) in our context, because we deal with the ring of holomorphic functions, 
but it is not usually considered in the CFT-context, where one prefers to work with the 
above rings of formal series. A related problem with our form of the consistency condition in 
the CFT-context is that the condition < |x — ?/| < \y\ < \x\, which is essential for our form 
of the consistency condition, does not even make sense in the ring of formal series considered 
in the CFT context. 

However, in the CFT context, there is another way to formulate a consistency condition 
which bypasses any convergence considerations, and which leads to a relation with a super- 
ficially similar appearance. This condition can be stated in various equivalent ways, but a 
particularly transparent formulation is that for each a,b E V there is an G N such that 

{x - y)^[Y{a, x)Y{b, y) - Y{b, y)Y{a, x)] = . (6.115) 

^ To be precise, our vertex operators should be written y(a,a;,a;), because they are not holomorphic in 
X unlike in the CFT context. 



43 



This equation now makes sense in the CFT context as an equation in the ring End(V^) ® 
C[[x,x~^,y,y~^]], but it no longer makes sense in our context: The first term would be 
defined in the domain |x| > \y\, whereas the second term would be defined in the disjoint 
domain \y\ > \x\. Thus, the situation with regards to our consistency condition fll.Sp and 
with fl6.115p is exactly opposite: Our condition (11 .Sp cannot even be formulated in the 
standard vertex algebra approach in CFT's, while eq. fl6.115p cannot even be formulated in 
our approach. While one can derive from eq. (16.1151) an equation that does make sense in 
the CFT-context, and that is similar in appearance to our consistency condition (11. 5p . 

{x - yf Y{a, X - y)Y{b, -y)c = (x - y)^ Y{Y{a, x)b, -y)c , (6.116) 

see [T^ Lemma 4.6], one sees that there are notable differences between this and our consis- 
tency conditioij^. 

The equation (I6.115p . which is often called "locality" in the CFT- vertex operator algebra 
context, may also be restated as saying that the commutator [Y{a, x), Y{b, y)] is a finite sum 
of derivatives of "formal delta distributions"!^ in the sense that 

[Yia,x),Y{b,y)] = 5^[9"5](a; - y)F(c„, , a,b,c^eV. (6.117) 

n=0 

Especially when written in the last form one can understand that the "locality condition" 
is intimately related to the fact that, in 2-dimensional CFT's, the commutator of local 
operators are supported on the lightcone (in a Minkowski formulation of the theory) — or 
equivalently — that the singularities of the OPE coefficients are of the form [x — ?/)~", where 
n G N (in a Euclidean formulation of the theory). By contrast, the singularities in a generic 
QFT do not have to be of this form, and may e.g. contain log's, as exemplified by the 
perturbative constructions of this paper. Furthermore, in a generic QFT, it is no longer true 
that the terms in the OPE are smooth apart from a finite number of terms, as even terms 
in the OPE 

Oaix)O,i0) = J2c:,ix) a(0) (6.118) 

c 

corresponding to fields Oc of arbitrarily large dimension still typically contain log's and are 
thus not strictly smooth. 

For these reasons, we are somewhat pessimistic that a fruitful definition of vertex algebra 
of the same purely algebraic flavor as described e.g. in [19] can also be found for the generic 

^In the approach [TTJ Sec. 3.4] to CFT's (see also references therein) the consistency condition is the 
same as ours. However, he does not impose the restriction < |a; — y| < |y| < for the vahdity of the 
equation. This feature is accidental in CFT and essentially due to the fact that three points 0,x,y € C 
can be brought into an arbitrary relative position by a conformal transformation, but not by a Euclidean 
transformation. Hence, also the underlying reasoning leading to the consistency condition as given in this 
reference is specific to CFT's and cannot be repeated in more generic situations. 

^°The formal delta distributions are defined as the formal series 6{x — y) = Xlnez a;"y^"^^. 
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(e.g., perturbative) QFT models. However, we think that this direction deserves further 
study, and some ideas in this direction have been put forward by [4J. 

Acknowledgements: We would like to thank R. Brunetti, J. Holland, N. Nikolov, and 
I. Runkel for discussions. We would also like to thank R. Brunetti for pointing out to us 
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A Spherical harmonics and Legendre functions in D 
dimensions 

Polynomials h{x), x G MP which are solutions to the Laplace equation Ah{x) = are called 
"harmonic polynomials". Since the Laplace operator A commutes with dilations x i— > tx, 
it follows that any harmonic polynomial can be decomposed into a sum of homogeneous 
harmonic polynomials. The harmonic polynomials satisfying h{tx) = t^h{x),l G N span a 
vector subspace of dimension N{l,D) in C[x], where N{0,D) = 1 and 

^,,B,.Pi±il_M_^ t-'>0. (A.n9) 

This can be seen for example by noting that the degree / harmonic polynomials h{x) are 
in one-to-one correspondence with totally symmetric, traceless tensors of rank / on M^: If 
C/xi...Ati are the components of such a tensor, then h{x) = X] c^ti...^j,a;^i ■ ■ - a;^, is a harmonic 
polynomial of degree /, and vice versa. The spherical harmonics in D dimensions are by 
definition the restrictions of the harmonic polynomials to S^^^. 

In the main text, we consider the a basis hi^mix), m G {1, . . . , A^(Z^, /)} of degree / 
harmonic polynomials for each / G N. The members of this basis are chosen to satisfy the or- 
thogonality condition eq. (13.281) . In fact, the hi^rn form an orthonormal basis of L'^{S^~^, dQ) 
when restricted to the sphere. It follows immediately from the fact that the hi^m{x) are har- 
monic polynomials that their restrictions hi^m{x) to the sphere are eigenf unctions of the 
Laplacian A on S^~^ with eigenvalue — /(/ + D — 2). 

For our calculations in appendix Oj we need to know in more detail the relation of the 
harmonic polynomials hi^m to the traceless symmetric tensors of rank / described above. To 
state the relevant facts, we use the familiar multi-index notation, a = (ai, . . . ,00) G N^, 
with 

a;" = J]a;°-, 9, = J] 9^-, a! = J]a^! etc., (A.120) 
M M M 



and we write 



hl,m{x) = J2km;aX'' . (A.121) 
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Combining eq. fl3.28l) with theorem 5.14 of [2] we get 

^ — \ - C(\ 

2_^'^l,m;a'tl' ,m';a-r = 5l,l'8m,m'- (A. 122) 

o. 

with ki = 2'r(/ + D/2)/T{D/2). This can also easily be proved starting from eq. (15.1091) . 
The decomposition of a harmonic function / regular at the origin into harmonic polynomials 
reads 

/(^) = Y.{j^^ dn{x)f{x)hi,^{x)\ hi^mix) . (A.123) 

l,m ^ ^ ^ 

With d"-x^\x=Q = Sa/sal we have 

hl,mid)hi'^rn' ix)\^=0 = ^ti^rn^aU' ,m';a(y\ = Su'Smm'kl (A.124) 



and thus eq. (IA.123P reads 



fix) = J2 {hiUd)f{0)) hUx). (A.125) 

l,m 

We also cite theorem 5.20 of [2], which states that for a harmonic homogeneous polynomial 
p of degree /, 

p{d) g{r) = qi r^~^p{x/r^) (A.126) 



where r = Ixl 



and 



( \- J ^^^^ for D = 2 
^^^^ ~ Mnr for D > 2, ( • ) 



. 2'-ir / for = 2 

qi = ■{ , ^ ' A.128 
^ 2'r(/ + D/2- l)/r(D/2- 1) forD>2. ^ ^ 



The Legendre polynomials in D > 2 dimensions are defined as the following invariants 
under SO{Dy. 

^(''^)_ 21 + D-2 

V hmi£)h,,miy) = P(x ■ y, I, D) . (A.129) 



m=l 



By construction, the Legendre polynomials P{z, I, D) are polynomials of degree / G N. They 
are often also called "Gegenbauer polynomials" and are denoted alternatively (z), 
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with notable differences in the normahzation convention throughout the hterature. A gen- 
erating function is 



^ . = y P(z, I, D)h' . (A.130) 

D-2 \y/l-2hz + hy t"o 

This formula holds for D > 3. For D = 2, the left side is to be replaced by 
— In \/l — 2hz + h'^. A generalization of this formula needed in the main text is pro- 
vided in theorem O The Legendre polynomials have the symmetry property P{z,l,D) = 
{—iyP{—z,l,D), and satisfy the normalization condition 

For complex values of the index z/ G C (or D), one can define an analytic continuation by 
means of the Gauss hypergeometric function 

Fiz, u, D) = u + D-2,D/2- 1/2, ^] . (A.132) 



T{u + 1)T{D-1 

The Gauss hypergeometric function is given by the convergent expansion 



n=Q 



for |x| < 1. Note that the above formula has a slight anomaly in D = 2 dimensions. Here, 
it gives P(cosa, z/, 2) = cos(z/a)/2z/ in D = 2 dimensions for u ^ 0, and this evidently does 
not have a limit as u ^ 0. On the other hand, the generating formula definition gives 
P(cosa,0,2) = 1. The differential equation satisfied by the Legendre functions is 

(1 - z^)y" -{D- l)zy' + u{u + D - 2)y = Q . (A.134) 



B Theorems for Legendre functions 

In the main text, we use certain identities for Legendre functions in D dimensions that we 
were not able to find in the literature, and which we therefore prove here: 

Theorem 1. {Generalized Dougall's formula) For any z/ G C \ Z and —l<z< +1 and 
-D > 3, we have the identity 



(B.135) 
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Proof: For D = 3, a. proof of the theorem can be given via a contour integral argument, 
see . We here give the following conceptually somewhat more transparent proof, valid for 
arbitrary D > 2. Let A be the Laplacian on the sphere S^~^. This is an elliptic, second 
order partial differential operator on a compact manifold with analytic coefficients. Using 
standard results on the functional calculus of such operators, we can form the resolvent 
operator = [A + v{v + D — 2)]^^ for any u such that z/(z/ + D — 2) is not an eigenvalue, 
i.e. z/ ^ Z. Let Ri,{x, y) be the kernel of Ry, which using general results on the Laplacian on 
compact Riemannian manifolds is known to be an analytic function on S^~^ x 5''^"^ apart 
from coincident points. Near coincident points, one has Ry ~ y)]"*^'^"^-'/^ for D > 3 
and Ry ~ \n.d{x,y) for D = 3, where d{x,y) = arccos(£ ■ y) is the geodesic distance on the 
sphere. A representation of Ry in terms of eigenfunctions of the Laplacian is 

oo NiD,l) - 

Ky[x,y) - 2^ 2^ ^^^^D-2)-l{l + D-2) 

1=0 m=l ^ J \ < I 

^ (2l + D-2)P{x-y,l,D) 

= "°l; „(„ + c_2)-i(; + c-2) - 

In the second line we have used the definition of the Legendre polynomials. Hence we see 
that the kernel Ry is, up to a constant, precisely equal to the left side of the Dougall formula. 

By definition, the kernel obeys [A + z/(z/ + D — 2)]Ry = S in the sense of distribu- 
tions. However, since Ry is evidently invariant under S'0(D)-transformations, we may write 
Ry{x,y) = y{z) for some analytic function of z = x ■ y when 2; 7^ 1. As a consequence of 
the differential equation satisfied by Ry, it can easily be seen that y satisfies the differential 
equation for the Legendre function of dimension D and degree u, see eq. (IA.134I) . Hence we 
have 

y{z) = A F{z, iy,D) + B F{-z, z/, D) (B.137) 

for some A,B ^ C as P{z,iy, D),P{—z,h', D) span the solution space of eq. (lA.134p . Fur- 
thermore, P{—z, u, D) is singular at 2; = 1 and regular at 2; = — 1 (see e.g. [8]), as is y{z). By 
contrast, P(z, D) is singular at z = —1 and regular at 2; = 1. Thus, we must have A = 
in eq. flB.137p . 

In order to determine the constant B, we evaluate y{—l). Using P{—l,l,D) = 
(-1)'P(1,/,D) and the formula flA.1311) . we find using various summation identities for the 
Gamma- function: 

J_ 1 f^. y/ 1 1 \ {l + D-3)\ 

an^^ ^ ^ \u - 1 u + l + D-2) l\ 

r(-z.)r(z. + D-2) _ T{v + D-2) TT 

r(D-i) r(z/ + i)r(D- i)sin7ri/ ■ ^ ■ ' 

Comparing this with the normalization of P{—z,v,D) at z = —1, we get the statement of 
the theorem. □ 
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The next theorem is a generahzation of formula (lA.130p . 

Theorem 2. {Shifted generating functional formula) For any 5 G C \ Z and —l<z< 
+1, \ h\ < 1 and even D > 4, we have the identity 

oo f 

J2 Pi^, l + 6,D) = i26Dy' AoUz + i sqrtl - 2F1 (Sd, 1; 6d + I; h{z + i sqrtl - z^)^ 



1=0 



+ {z-i sqrtl - z'^Y^ 2F1 {6d, 1; 1 + 5d, h{z - i sqrtl - z^) j | 

where 5d and the differential operator Ad are given by 

1 f d 



(B.139) 



T{D/2)[^2hdz) 
5d = S + {D-2)/2. (B.140) 

For D = 2 the differential operator is missing. For odd D > 5, we have the formula 

°° f -1 

h^F{z,l + 6,D) = Apl ^ -.X (B.141) 



1=0 



(-&.&. 1. 1; i^, H- (1 1 + 1. 2; i^. ) } 

where this time 

^ 2r(D/2) \2hdzJ 

5o = 6 + {D-3)/2. (B.142) 

For D = 3 the differential operator is missing. Fi is the two-variable generalization of the 
hypergeometric function defined by 



F,{a, b, c; d;x,y) = J2J2 ^ ,7"^ ^ , , , (B-143) 

t^atr'a {d)m+nm\n\ 



and we have defined t± = h ^ (l ± \'l + h?^—2hz) , and r(t) by eq. ( IB. 1551) . 



Remark: There is an apparent asymmetry in the formulas for even and odd D. One is 
tempted to believe that both formulas given for the shifted generating function are valid for 



49 



all D (when appropriately interpreted), but we have not been able to show this. The first 
formula may further be rewritten noting the standard formula 

f{5,x) = (25)-SFi(5,l;5 + l;l-a;) = ^^[V'(n + 1)- 7/^(5 + n)-lna;]a;" . (B.144) 

In the main text, we use the theorem to calculate the sums eq. fl4.7ip ioi D = 2 and eq. (14.931) 
for D = 3, where we also make use of the relation 

F{z, u, D) = (-l)^~2p(2, -u - D + 2,D) . (B.145) 

We then apply the recurrence relation (IB. 1461) to obtain (7d(5, cos/5, t) as in eq. (I4.94p (for 
even D). 

Proof for even D: For D = 2, the proof of the theorem follows immediately from the 
definition of the hypergeometric function. The alternate form ( IB. 1441) is obtained from 
transformation formula 15.3.10 of [T] for hypergeometric functions. For D even and D > 2, 
we prove the formula using the recurrence identity 

-^F{z,u,D)=DF{z,iy-l,D + 2). (B.146) 

Proof for odd D : For odd D, we proceed using the same recurrence identity, but in order to 
be able to do so, we have to evaluate Y1u=(d-?,)/2 ^ + 3), and this requires some extra 
work. We start with the Schlaefli integral formula for Legendre functions 



P(z, I/, 3) = — / ^dt . (B.147) 

To make (t^ — iy2^'^{t — z)^^~^ single- valued, we have to introduce two cuts in the complex 
plane, and we follow [26] choosing these cuts as the half-line 71 = (—00, —1) and a curve 72 
joining the points t = 1 and t = z, parametrized by 

(1< 71 < 00) . (B.148) 

The contour C"*" encircles 72 counterclockwise, see fig. [HI In [22], this particular representation 
of Legendre functions made it possible to determine derivatives of Legendre functions with 
respect to their degree u. 

We set 6d = 6 + {D - 3)/2 and have 

OD OD 

J2 h'P{z,l + 5D,3) = J2h'P{z,l + 5D,3) 

Z=(D-3)/2 /=0 

1 f . fe-iV^ 1 



dt 



2i+'5c7ri Jc+ \t- z J t- z- hif^ - 1) /2 



2^+^1^711 Jc+ 't- z-h{t'^ -l)/2 
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Figure 8: Cuts chosen to make (t^ — — z) ^ well-defined. 

where in the first equation, we have interchanged the order of summation and integration, 
and in the second we have set = (t^ — — z)~^° (defined with the same cuts as 

above). Our aim will be to express x as a sum of two functions XiyX2 that have cuts on 
7i and 72 respectively, and are analytic elsewhere. Then we will be able to carry out the 
integration in eq. (IB. 1491) . Assuming < 0, we apply formula 3.1.11 of [T7] . 



,w^-.-./d.*(s„.,.))(.-r'. (B.160) 

We introduce a function pi that equals 1 in a small neighborhood of the cut {—00, 1) and 
outside a slightly bigger neighborhood. Also, we introduce p2, smooth, equal to 1 in 
a small neighborhood of the cut flB.148p and in a slightly bigger neighborhood, so that 



supppi n suppp2 = 0- As X is analytic away from the cuts (i.e. 9x = on C \ (71 U 72)), we 
can modify eq. ( IB. 1501) in the following way, 

X{t) = Xi{t)+X2{t), 

Xi{t) = -n-' [ dwdw{dx{w,w))iw-t)-'p,iw,w), (1 = 1,2). (B.151) 



Now by theorem 3.1.12 of [IT] 



1 

Xi{t) = — da;(x(x + iO)-x(a:-iO))(a:-t)-i 
sm 0£)7r / X — 1^ 



^ ■' -00 \ ^ ^ / D 



dx\- {x-ty^. (B.152) 



We substitute m := 2/(1 — x) and obtain 



sin^DTT . /"^ A \-z\^'' \-t 



XI (t) = ^2'- j^u-'-'-{\-u%[\-^u^ (^^-^^ 

sm ^DT^_^^uB^^ + 1^ 1. (1 _ ^)/2, (1 _ t)/2) (B.153) 
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where in the second equation, we have used formula 3.211 of \TT\; Fi is the hypergeometric 
function of two variables defined above, and B is the Beta function B{x, y) = r{x)r{y) /r(x+ 
y). We next determine X2, 

X2it) = -TT^^ / dwdw{dg){w - ty^p2{w,w) 



-TT ^ j drdr^^ (^^dfx{w{T))^ (w(r)-t) V2(w(r), «;(r)) 

-— / dx {x{x + lO) - x{x - lO)) t 

2vri {z-xf \ z-x 

sm o/)7r / I X ~\\ 1 — z 1 — zx ^ 



I dxl ^ 1] (B.154) 

y_oo \z-xjjj{z-xy\z-x J 

where in the second equation we have performed a change of coordinates, 

t(w) = - —, (B.155) 

z — w 

with inverse w{t) = (1 — tz)/{z — r). This coordinate transformation maps 71 on 72 and 
vice versa. Again substituting u = 2/(1 — a;), we get 

sin5z3vrl-z2 (\-z \-^-^° ( \ - Tit) 

sin5fl7rl^^2_i+5^ 

TT t — Z 

X 5(l + 5z5,l-M^i(l-^D,l + 5D,l,2;(l-z)/2,(l-r(t))/2) (B.156) 
In eq. ( IB. 1491) . we replace 



t-z-h{t'^-l)/2 hit - t+){t - t_) 

t± = h-^ (l±Vl + h^- 2hz^ . (B.157) 

For — 1<-2<1, t+G7i and t^ G 72, so t- lies inside C"*" and t~^ outside. We can now 
calculate the contribution of to the integral ( lB.149p . 

y^£/''^'"' i-.-/,y-i)/2 - <2-i^«) 

which is now a simple residue integral. We obtain 

^ i?(l + 5z5,-Mi^i(-^D,5D,l,l;(l-^)/2,(l-t-)/2) . (B.159) 



TT VTTh^+2hz 
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In a similar manner, we calculate the contribution of X2, 

X / drf / ~ ""^ ^1(1 -6d,1 + Sd, 1, 2; (1 - z)/2, (1 - r)/2)- ^ 



— rj^ r(rj — z 



/ 2\ Z-T Z-T \ 

^ y—h) {t_-z){r-r{t_){U-z){r-r{U))^ ^^"^'^^ 

where we used the coordinate transformation (lB.155p . The contour D~ is the image of 
under this transformation. D~ encircles 72 and runs clockwise, not crossing the cuts 71,72- 
We can deform D~ into C~ by which we mean the contour with negative orientation. 
r(t+) is on the inside of C~, T(t_) on the outside. Thus the residue integral (IB. 1601) is 

sin tt5d 1 1 — z^ 



vr VI + - 2hz 2{z - t+) 
X B{1 + 6d, I -6d)F,{1-6d, I + 6d, 1,2; {1~z)/2,{1-t{U))/2). (B.161) 



Putting together eqs. (IB. 1591) and (IB. 1611) . we obtain 

sin ^^iTT 1 



1=0 



n Vl + + 2hz 
X ( 5(1 + 6d, -Sd)F, {-6d, Sd, 1, 1; (1 - z)/2, (1 - t_)/2) 



X Fi {1 -6d, I + Sd, 1,2; {1 - z) /2, {1 - t{U)) /2) ] . (B.162) 



Above, we have used < 0. However both sides of eq. (1B.162P are analytic in 6d through- 



out the complex plane for < 1 (with possible exceptions for (5/) G Z). This means that 
they are identical and eq. (IB. 1621) must hold for 9^5 > as well. Eq. ( IB. 1411) follows from the 
first line of eq. ( IB. 1491) . the recurrence formula (IB. 1461) . standard identities for the gamma 
function and the relation [z — t^){z — t^) = z^ — 1. □ 
We finally mention another representation of the Legendre functions used in the main 
text: 

Theorem 3. For z/gC\Z, |2;|<lwe have the formula 

V(..n\ sinvrz. 2-(^+^)/^ ^ o ^n r(-^/2 + n/2)r(z./2 + n/2 + D/2 - 1) 
nz,^.D) = —-^^^{-2z) . (B.163) 
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Proof: We prove this first for D = 2. Let z = cos a. We liave tlie identities 

cosi^a = 2^i(-i^/2,z//2; l/2;sin2a) (B.164) 

= -, ,^ , , --2Fi(-z//2,z//2; l/2;cos2a) 

r(-z//2 + l/2)r(z//2 + l/2) ^ /WW. ; 

-2 cos a y^_^/2)Tiu/2) ("''/^ + ^2, + 1/2; 3/2; cos^ a) , 

wliere in tlie second line we liave used a standard transformation formula for hypergeometric 
functions. We now use P(cosa, i/, 2) = cos(z/a)/2i/, and we expand the hypergeometric 
series in the second and third line, using the doubling identity of the Gamma function, 
\/7rr(2x) = 2^^~^r(x)r(x + 1/2), in various ways. Then we obtain the statement of the 
theorem for D = 2. The case D = 3 is covered by formula 8.1.4 of pQ, together with the use 
of the doubling identity as above. 

For general D G N, we use the recurrence formula (IB. 1461) . combined with a standard 
formula for the derivatives of the hypergeometric function. This then gives the formula for 
all even D starting from D = 2 and all odd D starting from D = 3. □ 



C The free field vertex operators 

The non-vanishing partial derivatives of the basic field ip in the theory defined by the 
Schwinger functions eq. (I3.30p are 

^''"^ = Cr\m;ad''^ (C.165) 

where q will be chosen later. The composite fields are labeled by multiindices 

l,m 

We use Latin letters for the multiindices denoting composite fields and Greek letters for 
multiindices when dealing with polynomials in x or d. The basic field if is harmonic by the 
field equation f l3.27p . so we may use eq. (1A.125I) . 

<^(x) = V^/i,,™(a;)<^''™(0) (C.167) 

1,171 

This has to be understood as an equation for insertions into Schwinger functions. Now the 
OPE of with a field Oa can easily be deduced from eq. fl3.30p and the definition of 99''"^, 
eq. flC.1651) : 
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where r = |x| as always; see eq. (1A.127P for the definition of the "Euchdean propagator" g. 
Again, this has to be understood as an equation for insertions: 

{ip{x)OaiOMxi) . . . ipixs)) = ((RHS of eq. Kum ) ip{xi) . . . ip{xs)) . (C.169) 

If we write down the left hand side of eq. flC.169p as in eq. fl3.3ip . where the points x and 
are represented by vertices v and w respectively ( = 0), then the first term of 

the right hand side in eq. (1C.168I) can be understood as a Taylor expansion (in x) of those 
graphs without a line (vw), and the second term as the Taylor expansion of graphs including 
that line. In order for these Taylor expansions to converge, we must have . . . , > 
Using eqs. (1A.126P and (1C.168P we obtain 

^(x)a(o) = V"''™ + ^ V ^^'^^^^ (0) 



Lm 



+ 5^v^^/i^,™(x)r-2'-^+2(^,_,^^^,(0) (C.170) 

where by e(;^m), we mean the multiindex defined by (e(/^m))z',m' = Si,i'Sm,m', and we define 
(a — := for ai^m = 0. To obtain a symmetric form of the OPE, we choose 

^— i2^T(l)^/T/2 forL) = 2 

ci = J^i= { ^ ^ ' ^ C.171 

[2'r(/ + D/2 - l)yj2{l + D/2 - 1)/{D - 2) for D > 2 . ^ ' 

We introduce the abstract vector space V spanned by the field labels a and creation and 
annihilation operators h^^^hi^rn on V by 



b^„,a = \/ai^rn + l(a + e(i,„)) 

b/,ma = ^ai.m{a - e(i,„)) . (C.172) 

We also introduce the vertex operator lo(v^) x) that corresponds to a multiplication of an 
insertion with the free field ^p{x). We rewrite the left-hand side of eq. ( 1C.170P in this notation 
by 

Fo(v',a;)a (C.173) 
and we can read off the right hand side of eq. (1C.170P that 

Fo(V^, x)=KoY. J— - {KMHrn + Km{x)r-''-''^^h,^n.) (C.174) 



with Kd = I ioi D = 2, Kd = for D > 2 and uj{D, I) = 21 + D - 2. 

That the vertex operators for composite fields Y^^a, x) are given by eq. fl3.35p is evident from 

the remarks on how to determine the Schwinger functions of composite operators from those 
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of the basic field below eq. fl3.30p . 

In the main text, we also used Wick's theorem, which we cite for completeness in the following 

form: Let bf^b^^, denote the "contraction" between two creation and/or annihilation 
operators, defined by 

Writing b* for b*^, and b* for b*^,, we have the combinatoric identity 

min(n,p) 

: bf ...b# : : b# ..b# : = ^ E bfb|-bTb| 

s=0 i\<...<is 

X :bf...b#b#...b#:(,„...,^,„...,,^) (C.175) 

where double dots mean normal ordering and the subscript (ii, z^; ji, j^) means that 
the respective creation or annihilation operators have been removed from the normal 
ordered product. 
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